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ABSTRACT 

We measure star formation rates (SFRs) of « 50,000 optically-selected galaxies in the local universe 
(z « 0.1), spanning a range from gas-rich dwarfs to massive ellipticals. We obtain dust-corrected SFRs by 
fitting the GALEX (ultraviolet) and SDSS (optical) photometry to a library of population synthesis models that 
include dust attenuation. For star-forming galaxies, our UV-based SFRs compare remarkably well with those 
derived from SDSS-measured emission lines (primarily Ha). Systematic deviations from perfect agreement 
between these two methods is shown to be due to differences in the dust attenuation estimates. In contrast to 
measurements based on Ha, UV provides reliable SFRs for galaxies with weak or no Ha emission, and where 
Ha is contaminated with an emission from an AGN (1/2 of the sample). We use full-SED SFRs to calibrate 
a simple prescription that uses GALEX far-UV magnitude and the UV slope to produce good dust-corrected 
SFRs for normal star-forming galaxies. The specific SFR (SFR normalized by stellar mass) is considered as 
a function of stellar mass for (1) star-forming galaxies with no AGN, (2) those hosting an AGN, and for (3) 
galaxies without Ha emission (the latter two groups forming the bulk of the optical red sequence). We find 
that the three have distinct star formation histories, with AGN lying intermediate between the star-forming 
and the quiescent galaxies. Normal star forming galaxies (without an AGN) lie on a relatively narrow linear 
sequence. Remarkably, galaxies hosting a strong AGN appear to represent the massive continuation of this 
sequence. On the other hand, weak AGN, while also massive, have lower SFR, sometimes extending to the 
realm of quiescent galaxies. We propose an evolutionary sequence for massive galaxies that smoothly connects 
normal star-forming galaxies to quiescent (red sequence) galaxies via strong and weak AGN. We confirm that 
some galaxies with no Ha emission show signs of star formation in the UV. We derive a UV-based cosmic star 
formation density at z = 0.1 with significantly smaller total error than previous measurements. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — ultraviolet: galaxies — surveys — 
galaxies: active 



1. INTRODUCTION 

Current studies of galaxies are characterized by two ma- 
jor features: availability of large samples of objects (galaxy 
surveys), and the utilization of the multiwavelength approach. 
Such studies extend from the nearby galaxies to those close 
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to the epoch of the formation of the first galaxies. A rate 
at which a galaxy forms stars is one of the more important 
properties in studying galaxy evolution. The multiwavelength 
approach allows us to employ a suite of star formation (SF) 
indicators — from X-rays to radio wavelengths. Major effort 
in recent years was made to provide reliable calibrations for 
different SF estimators and to understand their differences, 
advantages and limitations. The most straightforward evalua- 
tion is achieved by comparing two or more star formation in- 
dicators for the same set of objects. While new galaxy surveys 
provide large statistical samples with which one can attempt 
such studies, the sample of galaxies for which more than a 
single star formation indicator can be applied is not necessar- 
ily large. Moreover, various SF indicators are often applied 
for different redshift regimes. Global cosmic star formation 
history is therefore the result of studies that employ different 
SF indicators. 

Among the most frequently used star formation indica- 
tors are the UV continuum (usually at A < 2000 A), nebu- 
lar recombination lines (primarily Ha, but also [Oil]), far- 
IR d ust emission, and the synchrotron radio continuum at 2 1 
cm (lKennicuttlll983h . A compr ehensive review of most of 
these methods was presented in iKennicutl (11998 ). together 
with simple formulae for the conversion of the true flux into a 
star formation rate (SFR). Recently, the use of other SF indi- 
cators has been explored, such as the X -ray continuum (e.g., 
iDavid et al.lll992t iKilgard et alJl2002h . or the luminosity of 
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PAH features in the mid-IR dRoussel et all 1200 Star for- 
mation indicators are only as good as the assumptions that 
connect a certain observed luminosity to the actual star for- 
mation rate, and therefore all are sensitive to various system- 
atic uncertainties. One factor has proved more frustrating to 
account for than the others — the effect of dust obscuration on 
the UV and Ha star formation rates. In contrast to these two 
dust-sensitive indicators, bolometric IR luminosity and the ra- 
dio luminosity are often considered to be "true" SF indicators. 
While this is certainly an oversimplification, such perception 
is bolstered by the very tight correlation of the IR and the ra- 
dio luminosities, at le ast as observed for the normal galaxies 
in the local universe (Ide Jong et al.lll985h . However, even if 
they were perfect, IR or radio methods cannot make UV or 
Ha methods obsolete for many practical reasons. 

There have been a number of studies that compared the 
UV and Ha star formation rates — either as the observed lu- 
minosities, or by applying various schemes to correct for the 
dust attenuation.'^ These comparisons were of ten carried out 
with respect to some dust-free SF indicator. Hopki ns et al.l 
(|2001) used a sample of several hundred objects with Ha and 
UV measurements (with the U -band photometry serving as a 
proxy for the UV), and compared them to the far-IR measure- 
ments compiled bv lCram et al.l ( Il998h . They find that in order 
to reconcile Ha and far-IR luminosities one needs to apply a 
dust attenuation that is not fixed, but rather depends on the 
SFR itself. Actually, this dependence is a consequence of the 
relationship betwe en the dust attenua t ion an d the far-IR lu- 
minos ity found by lWang & HeckmanI (Il996l) . iHopkins et alJ 
(1200 Ih use an extinction curve to extrapolate the dust attenu- 
ations obtained for Ha into the UV regime, but find that such 
correction fails to bring UV luminosities to agree with the far- 
IR, i.e., the simple application of the same attenuation mech- 
anism to b oth the Ha and the UV em ission does not appear to 
be correct. Bell & KennicuttI (1200 lb used actual UV observa- 
tions of 50 nearby galaxies in two ultraviolet bands obtained 
with the Ultraviolet Imaging Telescope (UIT), and compared 
them to Ha luminosities from the narrow -band imaging. They 
determine UV attenuati on for 13 galaxies using the correla- 
tion with the UV slope (ICalzetti et alJll99 4?). and measure Ha 
attenuation for 21 objects using the thermal radio continuum 
fluxes. They find that both can reach high levels even for nor- 
mal galaxies (~ 4 mag for Ha, and > 5 mag for far-UV), and 
lend suppo rt to previous notions that th e attenuation increases 
with SFR. Sullivan et all (l2000l 1200 lb confirmed that a bet- 
ter agreement be tween the UV obser vations from the FOCA 
balloon mission ( Mi Uiard et al]fT992l) . and the fiber Ha spec- 
tra is achieved when attenuation corrections are taken to be 
luminosity-dependent. They also suggest that one perhaps 
cannot use simple extinction curve scalings to convert Ha 
attenuations into UV attenuations. The breakthrough in re- 
solving this problem came with the i ntroduc tion of the two- 
component dust attenuation model of ICharlot & Fall (.2000) . 
This model was indeed motivated by the need to produce a 
consistent model for dust attenuation affecting Ha and UV 
continuum photons. It postulates the existence of short-lived 
(10 Myr) giant molecular clouds that affect photons produc- 
ing the Ha line. On the other hand, the attenuation of the 
UV continuum, having timescales longer than the Ufetime of 

In this paper we will use the term "attenuation", rather than the more 
commonly used term "extinction", thus emphasizing the complex processes 
of absorption and scattering in a galaxy, rather than the dimming of light 
along a single line of sight. 



giant molecular clouds, is predominantly produced by the dif- 
fuse ISM (after the molecular clouds have dispersed), at levels 
that are typically 3 times lower (for a given wavelength) than 
tho se in the molecular clo uds. In this paper we will be using 
the ICharlot & Falll (l2000h model, thus testing it for the first 
time on a large scale. 

In addition to systematic trends, we should mention that 
some previous studies were finding that the measurement 
errors are smaller than the observed scat ter. ISuUivan et alJ 
(2000V2001b: IIglesias-Paramo et all (12004 ') offer the explana- 
tion for this scatter as arising from the differing timescales 
over which Ha and UV SFRs are sensitive, so that UV could 
in some cases (especially in low-mass galaxies) detect a short 
starburst that is no longer observable in Ha. 

Sometimes, the systematics are present in the observations 
themselves. Rosa-Gonzalez et al. (2002) showed that even the 
relatively reliable determination of Ha attenuation can be af- 
fected by the systematics when not correcting for the under- 
lying stellar absorption in Balmer emission lines. Systematics 
can also arise wi th galaxy samples selected at different wave- 
lengths dBuat et a l. 2002). 

The main obstacle in obtaining UV measurements for a 
large number of galaxies in the local universe is that they 
need to be made outside of Earth's atmosphere. For this rea- 
son in 2003 NASA launched the Galaxy Evolution Explorer 
(GALEX, Martin et al. 2005). GALEX is currently conduct- 
ing the first ever survey of the UV sky. The imaging (in two 
UV bands) is executed in several modes — from a shallow all- 
sky survey, to the ultra-deep fields. In this paper we are us- 
ing measurements obtained in the medium-deep survey that 
is designed to image regions of the sky covered by the Sloan 
Digital Sky Survey (SDSS). Thus we obtain a large sample 
of galaxies with both UV and optical photometry, as well as 
spectroscopic redshifts. GALEX and SDSS data, and the re- 
sulting sample, are described in §^and[3] Dust-corrected 
star formation rates (and some other physical properties such 
as the stellar mass) are obtained by comparing the observed 
colors to stellar population synthesis models to which the dust 
attenuation has been applied (®. While we use full UV to 
optical SED, the SFRs are essentially constrained by the UV 
(Salimetal. 2005). The large size of our sample permits a 
robust statistical analysis. We com pare our UV-based SFR s 
to the results of the major study of Brinchmann et aP (l2004l) . 
who use SDSS spectra to derive Ha-based SFRs for ~ 10^ 
galaxies (©. They employ an intricate scheme to correct for 
the fiber aperture effects. Thus, our study also serves as a 
check on the reliability of their methodology. Finally, we dis- 
cuss star formation histories of different classes of galaxies 
(© and derive a UV-based estimate of the global star forma- 
tion density at z = 0.1 (®. 

2. DATA 

2.1. GALEX data 

Technical aspects of G ALEX telescope , detectors an d data 
products are presented in'Morrisse y et alJ (l2005Ll2007h . Here 
we give a summary of relevant information. GALEX surveys 
the sky in either the imaging or the grism spectroscopy mode. 
It simultaneously produces a far-UV (FUV) and a near-UV 
(NUV) image having a circular field of view of 1.25° diame- 
ter. FUV and NUV filters are broadband, with effective wave- 
lengths of 1528 A and 2271 A, respectively. We will denote 
the magnitudes measured in these photometric bands as FUV 
and NUV. A single field of view imaged by GALEX is called 
a tile. A tile can consist of one or more visits, i.e., individual 
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exposures. GALEX surveys the sky in several imaging modes, 
which differ in the exposure time per tile. In this study we 
use the Medium Imaging Survey (MIS), which is designed to 
maximize the coverage of the sky that is included in the SDSS. 
Typical exposure times in MIS are 1500 s, yielding limiting 
magnitudes of FUV = NUV = 22.7 mag (AB system is used 
throughout). GALEX cannot point in the vicinity of bright 
sources (usually stars), which inevitably produces some gaps 
in the coverage. 

GALEX data used in this paper come from the MIS por- 
tion of the Internal data release 1.1 (IRl.l), which is an ex- 
panded version of the first public GALEX release (GRl) of 
MIS (in both IRl.l and GRl the same pipeline, version 4, is 
used to reduce the data and produce source catalogs). Dataset 
consists of 705 GALEX tiles. Because of the anomaly with 
the FUV detector, 98 tiles lack FUV images. The 705 tiles 
cover 797 sq.deg. of the sky. Each GALEX field is restricted 
to 0.6° radius, close to the maximum field of view.''* Because 
of the overlap between GALEX tiles, the total unique area is 
741 sq.deg. Source catalog for each tile was produced using 
the SExtractor software (Bertin & Arnouts 1996). 

In our study we use FUV and NUV fluxes measured in Kron 
elliptical apertures. We recalculate flux errors because they 
were incorrect in th e pipeUne version 4 r eductions (thus also 
affecting GRl, cf. ITrammell et al.l 12007 ). Kron magnitudes 
should represent a good measure of a galaxy's total flux. If 
both the FUV and the NUV detections are present, we use 
the aperture (size, shape and position) derived from the NUV 
image to measure the FUV light. We go from NUV to FUV 
because we are generally more sensitive in the NUV. We add 
zero-point calibration errors of 0.052 and 0.026 magnitudes to 
FUV and NUV respectively to account for systematic inaccu- 
racies. The calibration errors were estimated by analyzing the 
repeat imaging of a calibration stellar object, and confirmed 
by a la rge number of repeat observations of the same field 
(iMorris sey et al. 2007). If a UV detection is present in only 
one band, we still measure the formal flux and its error in the 
other band (using the aperture defined in the detected band). 
In some cases, the optical source is not detected at all by 
GALEX, in which case we note the sky background at the po- 
sition of the SDSS source and compute the flux error. In this 
calculation we use galaxy sizes that have been derived (us- 
ing a calibration based on objects detected by GALEX) from 
SDSS Petrosian radii in r-band. Finally, in some cases the 
FUV cannot be used because the FUV image does not exist 
(15% of all galaxies). 

2.2. SDSS data 

In addition to the ultraviolet data from GALEX, we 
use o ptical d ata from SDSS Data Release 4 (DR4, 
lAdelman-McCarthv et al.1 l2006h . SDSS is providing five- 
band broadband photometry (ugriz bands), and the spec- 
troscopic follow-up of most galaxies with r < 17.77 (main 
galaxy spectroscopic survey. IStrauss et ani2002b . In addition 
to the official data products from the SDSS collaboration, 
we use value-added galaxy catalogs produced by MPA/JHU 
SDSS team. These catalogs include reprocessed SDSS spec- 
troscopic data, and some derived galaxy physical parameters 

'* Many other studies with GALiiX restrict analysis to 0.5 or 0.55° radius 
field of view. This is because artifacts are more common near the detector 
edge, and the PSF becomes distorted. Since in this study we match SDSS 
objects to GALEX sources, the chances of a match with an artifact are small. 
Also, while the P SF a t the edges is distorted, the total GALEX flux is not 
affected (see also j)4.3) . 



that are based on the spectroscopic data. 

SDSS catalog lists magnitudes measured in several differ- 
ent ways. We use MODELMAG magnitudes, which are the 
measurements of choice for relative fluxes (i.e., provide sta- 
ble colors), while still capturing most of the total light. We 
apply slight adjustments (-0.04 and 0.02 mag) to u and z 
magnitudes respectively to bring them closer to the actual 
AB system jAbazaiian et al.ll2004l) . We convert SDSS mag- 
ni tudes and errors into fl uxes using the transformations given 
in IScranton et al.l (l2005h . In addition to catalog photometric 
errors, we add 0.01 mag of uncorrected calibration error to 
each of the bands (see 34. 3t , and, specifically for the u band, 
we add an additional color-dependent error d ue to the red leak 
(o-«.RL = 0.0865(r-0^+0.0679(r-/), based onl Abazaiian et alj 
rare cases an SDSS magnitude would be missing 
and such band has to be excluded from the analysis. We fur- 
ther exclude from analysis individual SDSS magnitudes that 
are nominally fainter than ugriz = 25. These are invariably 
spurious, regardless of the listed photometry error. Finally, we 
exclude magnitudes with errors larger than 10 x the typical er- 
ror in that band (we find typical errors to be 0.07, 0.009, 0.007, 
0.007 and 0.017 mag for u, g, r, i, and z bands respectively). 
While large, we find these errors to be significantly underes- 
timated. To summarize, except in cases listed above when we 
exclude certain individual flux points, our input data consists 
of 7-band photometry and spectroscopic redshifts from SDSS. 

2.3. GALEX/SDSS matched catalog 

Of 741 sq.deg. of GALEX unique imaging, 645 sq.deg. 
overlaps with SDSS DR4 spectroscopic area, thus defining 
the solid angle of our sample.'^ We estimate this area by 
counting the SDSS galaxies that: (1) have spectra, (2) have 
a dereddened magnitude 14.5 < rpetro < 17.5 (the faint end 
is taken to be comfortably brighter than the spatially variable 
spectroscopic limit) and (3) lie at redshifts 0.005 < z < 0.22. 
We count such galaxies both in our survey, and in the full 
DR4, whose spectroscopic area is known. The ratio of the 
two counts gives us the size of our survey area. This estimate 
should be good to within 1 sq.deg. There are 67,883 objects 
from the SDSS DR4 spectroscopic survey (not restricted to 
the main galaxy survey) that lie in this area, and are spectro- 
scopically classified as galaxies. For each of these objects we 
search for a match in the GALEX source catalog (which al- 
ready combines FUV and NUV detections) within 4". Our 
analysis of SDSS point sources with a match in GALEX in- 
dicates that GALEX positions have a random error of 0."8 in 
either R.A. or declination (becoming somewhat larger at the 
edges of the field). In addition, there are overall tile-to-tile off- 
sets between GALEX and SDSS coordinate systems of several 
tenths of an arcsecond. In any case, astrometric uncertainties 
are significan tly smaller than our matching radius (see also 
ITrammell et al. 2007). A genuine match can be missed if the 
centroid of the optical light is offset by more than 4" com- 
pared to the centroid of UV light. We expect such cases to be 
quite rare, since at the mean redshift of the sample (z = 0. 104) 
this offset would translate into a 7 kpc separation. A problem 
with matching in general is that what is considered to be a 
single object in one catalog can be resolved into multiple ob- 
jects in another catalog, whether they are indeed separate ob- 
jects (blending), or actually belong to the same system (shred- 

Matching with the current SDSS data release DR5 (which uses the same 
pipeline as DR4) would not significantly increase the overlap, since current 
GAL£X pointings mostly follow the footprint of SDSS DRl and DR2. 
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ding). This problem is more pronounced when combining cat- 
alogs with different resolutions, as is the case here (4 to 5" for 
GALEX vs. 1"4 for SDSS). However, in our particular sam- 
ple this problem is not severe. SDSS galaxies with spectra are 
relatively bright objects, and if more than one GALEX object 
is found within the search radius our procedure was to take 
the one that is positionally the closest. Since we are dealing 
with high-latitude fields where the density of foreground stars 
is not that high, this simple matching procedure produces re- 
liable results. Besides, we do have mechanisms of identifying 
the majority of incorrect matches at the later stage, when we 
perform an S EP fi tting to the combined GALEX/SDSS pho- 
tometry (see ^4.3b . Since we combine photometry from two 
different catalogs there is a concern of photometric zero-point 
mismatch. We explore such a possibility in ^5.41 In our final 
matched catalog we eliminate duplicate GALEX observations 
(stemming from overlaps or repeat observations) by keeping 
those that lies closer to the center of GALEX field of view. 

3. THE SAMPLE 
3.1. Sample selection 

We initially define our sample by applying magnitude and 
redshift cuts to galaxies with SDSS spectroscopy (note that 
objects spectroscopically classified as QSOs are excluded). 
We require the dereddened Petrosian magnitude to be in the 
14.5 < rpetro < 17.77 range, where the faint limit is the nom- 
inal limit of the SDSS main galaxy spectroscopy survey (see 
also ^4.41) . while the bright limit is chosen to avoid objects 
with saturated SDSS photometry. We require redshifts to lie 
within 0.005 < z < 0.22 range. The lower limit is chosen to 
eliminate galaxies that may deviate from the Hubble flow, i.e., 
whose redshift distances could be unreliable. Redshift dis- 
tribution in SDSS main galaxy sample peaks sharply around 
z = 0.1 with few galaxies beyond our upper redshift limit. 
T hese magnitude and redshift limits are identical to those used 
in iBrinchmann et al.l (12004 ). to which we will be comparing 
many of our results. We will refer to the above redshift range 
as representing the "local universe". 

There are 49,346 galaxies that meet the conditions defin- 
ing our initial sample. Good quality SED fitting (see 34.2l i 
is obtained for 98% of them (48,295 galaxies). In the rest of 
the paper we will use only this latter sample and refer to it as 
"all" galaxies. Note that since we retain objects regardless of 
whether they were detected by GALEX (as long as they fall 
within the area covered by GALEX), our sample is only opti- 
cally (r-band) selected. 

Throughout the paper we will use fl,,, = 0.3, Ha = 0.7, Hq = 
70kms"' Mpc"' (i.e., hjo = 1) cosmology. 

3.2. Emission line diagnostics and sample classification 

This work focuses on the physical properties of galaxies, 
such as their star formation rates (SFRs) and stellar masses. 
Before we start analyzing the sample based on these derived 
properties, we would like to characterize it in terms of its ob- 
servable quantities. 

Throughout this paper we will rely on optical emission lines 
to classify galaxies in our samp le. This classification, based 
on the BPT diagram dBaldwin et al.lll98lb . plots the flux ra- 
tio of [OIII]A5007A and H/3 lines against the flux ratio of 
[NII]A6584A and Ha li nes. In this paper we fully adopt 
the BPT classification of Brin chmann et al.l (2004) (hereafter 
B04), which is based on emission-line strengths corrected for 
the underlying stellar absorption (see their Fig. 1). The posi- 



tion of a galaxy in the BPT diagram indicates the nature of its 
ionizing source. The majority of galaxies in the diagram fall 
within one of the two branches. One branch lies mostly above 
the maximum line ratios expected from star formation. These 
galaxies must have some ionizing source in addition to young 
stars. Most researchers attribute this emission to an active 
galactic nucleus (AGN). Specifically, the emission is associ- 
ated with a narrow-line AGN (a LINER or a Sey fert 2), also 
known as Type 2 AGN (iKauffmann et alJl2003d) . Note that 
any Type 1 AGN spectrum (QSO or Seyfert 1) would have 
broad lines, and would be classified by the SDSS pipeline 
as a QSO and thus not i ncluded in this sample. Following 
IKauffmann et alJ (l2003ch . B04 classify galaxies that lie on 
the lower portion of the AGN branch, but with line ratios al- 
lowed by star formation, as "composite" (star forming/AGN) 
galaxies. For galaxies that lie on the star forming branch we 
expect very little emi ssion line flux to come from an AGN 
(IKauffmann et alJl2003G) . B04 required a S/N ratio in all four 
lines to be > 3 in order to apply the BPT diagram classifi- 
cation. However, in some cases a secure AGN classification 
is possible even when only [NII]6584 and Ha have S/N> 3, 
i.e., in cases when this ratio is larger then the one allowed 
by SF B04 cafls this class a "low-S/N AGN (LINER)". Fol- 
lowing B04, we will study this class together with the AGN. 
B04 also introduce the category of "low-S/N SF" galaxies. 
These are the galaxies that do not fall in any of the previous 
categories because their lines have low S/N, yet they have an 
Ha detection with S/N> 2. While their lines are not strong 
enough for secure placement on the BPT diagram, if we nev- 
ertheless do so, we find that many object lie in the high-mass 
end part of the SF sequence, as well as in the lower portion of 
the AGN branch. Therefore, these objects represent a hetero- 
geneous class. Finally, there are galaxies without detectable 
lines, thus precluding the classification in the above scheme. 
B04 call this group "unclassifiable". We will call them "No 
Ha" class. We find that in ~ 2% of galaxies in this class the 
Ha non-detection is due to some artifact in the spectrum or 
line-extraction pipeline. We exclude these galaxies from this 
class (but not from the whole sample). The fraction of galax- 
ies in our sample belonging to different classes is as follows: 
Stai- forming (SF) - 27%, Low S/N SF - 19%, Composite 
(SF/AGN) - 8%, AGN - 12% and No Ha - 33%. 

3.3. UV-to-optical color-magnitude diagram 

Color-magnitude diagram (CMD) is a powerful tool in as- 
sessing the basic properties of a sample of galaxies. Histori- 
cally, the study of CMDs was preceded by the studies of a re- 
lationship between color and morphology. Basic morpholog- 
ical segregation into disk-like spiral galaxies and sp heroidal 
elliptical galaxies was established bylHubblel (Il926l). After- 
ward (e.g., Hubble 1936) it was realized that spiral (late type) 
galaxies have bluer colors than the ellipticals (early types). 
Optical CMDs were con structed for cluster early type galax- 
ies (e.g., ISandagdll972l) . where they were found to form a 
narrow sequence (the so called red sequence). Field spi- 
ral galaxies also displaye d the color-magnitude relationship 
dChester & Robertsl ll964). albeit with larger scatter The bi- 
modal nature of the field galaxy CMD became much more 
apparent with recen t large scale surveys, in particular SDSS 
dBaldrv et alj|2004b . However, unlike the traditional optical 
CMDs, a CMD in which the color is constructed from an ul- 
traviolet and an optical magnitude has a particular diagnostic 
power By virtue of contrasting the current (or recent) star for- 
mation as indicated by the UV light (modulo attenuation) to 
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Fig. 1. — UV to optical color-magnitude diagrams (CMDs). Upper left panel shows a greyscale scatter plot of all galaxies in our sample with a near-UV 
detection. (The shade of gray is directly coirelated to the number of points contained in a given "pixel".) Note the pronounced bimodality of the blue and the red 
sequences, and their large separation. Dashed line represents the completeness limit at the mean redshift of our sample. The remaining panels show CMDs of 
different galaxy classes (SF - star forming, Comp - star-forming/ AGN composite, AGN - Type 2 AGN, SF (low S/N) - star-forming with weak Ha, and No Ha), 
as determined from the position in the BPT diagram. The outer contour encompasses 90% of the entire sample is plotted for reference. While SF galaxies mostly 
lie in the blue sequence, and those with no Ha in the red, most galaxies in between the two sequences are AGN or AGN/SF composites. Absolute magnitude is 
given in z = r-band, and the color is K-corrected to z = rest-frame, as indicated by superscript 0. 

the total past star formation as indicated by the optical light, 
the UV to o ptical color is a good proxy of a galaxy's SF his- 
tory (e.g. JSalim et alj|2005b . On the other axis, one plots ab- 
solute optical magnitude, which is related (modulo variations 
in the optical mass-to-light ratio) to galaxy's current stellar 
mass. We present the observed CMD of our sample in Figure 
[1] In this and many subsequent figures, the individual data 
points have been converted into a greyscale density scatter 
plot, in such a way that the shade of gray is proportional to 
the number of objects occupying a given "pixel". Such repre- 
sentation is desirable when dealing with large samples where 
it is easy to saturate a traditional scatter plot. Unless noted 
otherwise, the figures are constructed from raw, unweighted 
counts. We will include volume corrections later, when ap- 
propriate. Both the color and the absolute magnitude have 
been K-corrected to z = rest-frame bands ( ^4.4b . The up- 
per left panel shows the CMD of all galaxies with an NUV 
detection (85% of the total sample; see Table [T] for a break- 
down of UV detection rates per class). A striking feature of 
a UV to optical CMD is the pronounced bimodality of blue 
and red galaxies. Blue galaxies form a well-defined sequence 
extending to faint luminosities. The red sequence is some- 
what more narrow than the blue sequence (note that this be- 
comes evident only after K corrections have been applied), 
and extends to intrinsically more luminous galaxies than the 
blue sequence. This is related to well-known fact that the 
most massive elliptic al galaxies are m ore massive than the 
most massive spirals (lHolmberglll965h . The two sequences 
are separated by ~ 3 mag. If the two sequences are modeled 
as gauss ians, one finds tha t there is an excess of galaxies in 
the gap (IWyder et alJl200'7h . This is not the case in classical 



All= 


48295 


86% 


0.38 


0.43 


SF 


12901 


99% 


0.20 


0.29 


SF (low S/N) 


9060 


93% 


0.30 


0.39 


Comp 


3966 


96% 


0.28 


0.40 


AGN 


5827 


90% 


0.41 


0.49 


No Ha 


16159 


68%, 


0.60 


0.54 



UV-based SFRs are averaged over 100 Myr.*" Either FUV or NUV detection. 
Of objects detected in FUV, 98% are also detected i n N UV-' Includes 382 
objects for which classification was not possible, see j|3.2l 



optical CMDs (iBaldrv et al.ll2004l) . We will ref er to the gap 
region and its population as the "green valley" dMartin et al.l 
12007b . A detailed quantitative study of the G ALEX UV to 
optical CMD is presented in iWvder et al.1 (l2007h . 

In the subsequent panels in Figure [T| we display the CMDs 
of various classes of galaxies as defined in ^3.21 Each panel 
has the greyscale normalized to the number of galaxies in the 
given class. For better reference with respect to the full sam- 
ple, in each subsequent panel we repeat the contour contain- 
ing 90% of all galaxies. Not surprisingly, the star-forming 
galaxies (SF) occupy the blue portion of the CMD. Part of the 
width of the sequence is due to the intrinsic dust attenuation. 
The CMD of composite galaxies (showing both signatures 
of SF and AGN) are shown in the upper right panel. Their 
NUV - r colors are offset to the red compared to those of the 
"pure" SF galaxies. Galaxies with narrow line AGN (lower 
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left panel) occu py the regions of the red sequence and o f 
the green valley (iMartin et alJl2007l: iKauffmann et alJl2007h . 
Moreover, most galaxies with intermediate colors are either 
AGN or composites. Star formation histories of AGN will be 
discussed more extensively in 37.31 Next, we have low-S/N 
SF galaxies (middle lower panel), which indeed mostly lie on 
the blue, star-forming sequence (preferentially its more lumi- 
nous part). However, there is a tail of red galaxies which are 
probably contaminated by AGN as observed in 33.21 Finally, 
we have galaxies with no detectable Ha (lower right panel). 
As expected, these galaxies form the bulk of the red sequence. 
However, there is a tail of galaxies of this class extending to- 
wards the blue colors. Since the spectral classification used 
here is nuclear (central 3") it is possible that some of these 
galaxies are classified as not having Ha because of the domi- 
nant bulge, while the relatively faint star forming disk is giv- 
ing a galaxy an overall blue color However, as we discuss in 
37.41 this is not the case for most of them. The CMDs pre- 
sented in Figure [Urequire a detection in NUV. The remaining 
15% of our sample with no NUV detection falls mostly in the 
red sequence (as evidenced from their u- r colors). These 
galaxies are too faint to be detected in the UV. 

4. OBTAINING GALAXY PROPERTIES WITH SED FITTING 
4.1. SEDs of model galaxies 

Spectral energy distribution (SED) fitting is becoming 
a widely used technique for de riving galax y prop erties. 
It was pioneered in w orks of ISearle et"aLl (11973 ) and 
iLarson & TinslevI (Il978l) . To first order it consists of com- 
paring the observed SED to a set of model or template SEDs, 
and searching for the best match. It is assumed that the phys- 
ical properties of model or template galaxies are known, and 
that one can use this knowledge to deduce the properties of an 
observed galaxy. 

In this study, we compare the observed SED to a large 
number of model SEDs constructed from Bruzual & Charloj 
(|2003) population synthesis code. To construct the models 
we first choose input model properties by randomly select- 
ing values from prior distributions defined in the following 
way. Formation of a galaxy is uniformly distributed between 
0.1 Gyr and the age of the universe. Global metallicity is 
uniformly distributed between 0.1 and 2 Zq. Star formation 
histories are not single stellar populations, but the combi- 
nation of an exponentially declining continuous star forma- 
tion, SFR(f) oc exp(-7f), with < 7 < 1 Gyr"' uniformly dis- 
tributed over that range, and the random starbursts superim- 
posed on the continuous SF. Bursts are parametrized to have a 
duration uniformly distributed in the 30-300 Myr range, with 
a strength such that the mass produced in a burst is between 
0.03 and 4 times the mass produced in the continuous SF over 
the present lifetime of a galaxy (distributed uniformly in log). 
Such parametrization conforms with most observational stud- 
ies. Finally, bursts are randomly produced so that the proba- 
bility that a given galaxy has undergone at least one burst over 
a 2 Gyr period is 50%. The above parameters define how the 
population synthesis for each model will be carried out. Next, 
each of the above model SEDs is subjected to dust attenua- 
tion parametrized acc ording to the two-component model of 
ICharlot & Falil ( l2000h . For the V-band optical depth we ran- 
domly chose a value from a distribution < tv < 6 (it peaks 
at Ti/ = 1 .2 and has most values in the < Ty < 2 range). The 
prior distribution for Ty is empirical, and comes from Balmer 
decrements in SDSS spectra. The choice of Ty prior distribu- 
tion will be discussed in 35.41 Fraction of the optical depth 



that affects stellar populations older than 10 Myr (i.e., most 
of the UV continuum flux) is given by factor /i, with values 
ranging from 0.1 to 1, peaking around /i = 0.3. Altogether, 
we produce 100,000 model spectra spanning a range of SF 
histories, metallicities and dust attenuations. Note that since 
we pick input parameters randomly, we do not call our set of 
models a grid, which would suggests a set of points with equal 
spacing in parameter space. 

Finally, we convolve the resulting model spectra with 
GALEX and SDSS bandpasses at five redshifts equally spaced 
in the [0.05, 0.25] interval, producing the libraries of model 
broad-band photometry. In each library we keep only mod- 
els that have an age smaller than the age of the universe at 
that redshift. This effectively reduces the number of models 
from 95,000 at z = 0.05 to 78,000 at z = 0.25. We also add 
effects of the interv ening intergalactic extinction according to 
iMadau et al.l (Il996h . Our final libraries list model photom- 
etry as well as a number of galaxy parameters associated to 
that model (such as the SFR averaged over several timescales, 
stellar masses, dust attenuation parameters, etc.) 

4.2. SED fitting 

In our study, the observed SED is constructed from GALEX 
and SDSS broadband photometric fluxes (the broad-band 
SED). For an observed galaxy at some redshift we select the 
model library with the closest redshift. We step through a li- 
brary one model at the time. Model flux points will have some 
arbitrary zero point, and in order to see how well an observed 
flux compares to the model flux we first need to find a factor 
a that minimizes the between the observed and the model 
points. In other words, for model /, we need to minimize the 
following expression: 



2 sr^ f Fohs.x-aiF^odi.x 
Xi=2^' — 



X 



\ o-(Foh!.,x) 



(1) 



Here, the summation is over X = (FUV,NUV,u,g,r,i,z), the 
observed flux points are Fo\,^, and their errors a{Fohs)- Flux 
points of model / are Fmod/ ■ By taking the derivative of Equa- 
tion[T]with respect to a, and equating it with zero, we find the 
scale factor a* that best matches the observed and the model 
flux points: 



X cr-(F„bs,x) 



(2) 



Because we include non-detections that have meaningless 
fluxes (small or even negative values), we use the observed 
fluxes without correcting for Galactic reddening, and instead 
apply the reddening to model fluxes, which are always pos- 
itive. This allows us to treat non-detections like any other 
flux point with a known flux error Reddening corrections for 
SDSS magnitudes come from SDSS d atabase, and those for 
GALEX from relations in lWyder et aP (l2007h (Their Eq. 2 for 
NUV). 

The goodness of the fit between an observed galaxy and 
model / is then obtained by inserting a* in place of a, in Equa- 
tion [U Note that only a single parameter is being fit — the 
scale factor a. We emphasize that the word "parameter" when 
fitting is discussed should not be confused with (in princi- 
ple) arbitrary number of galaxy parameters (i.e., properties) 
that correspond to each model SED. Therefore, since we have 
seven photometric points, there are six degrees of freedom. 
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For a given galaxy we evaluate xl for each model as out- 
lined above. Now, each of these xj values will determine 
the relative weight w, for the given model / (and therefore the 
weight of the galaxy parameters associated with that model) 
as w, = exp{-xf /2). We then build a probability distribution 
function (PDF) for every galaxy parameter of interest by com- 
pounding weights at the corresponding parameter value. We 
repeat the procedure for each of the 10^ models in the given 
redshift library. We normalize the final PDFs and note the pa- 
rameter values corresponding to the 2.5, 16, 50 (median), 84 
and 97.5 percentiles of the cumulative PDF. We use the av- 
erage of the PDF as our nominal estimate of the parameter 
value. Since in the general case the PDF is not symmetric, 
medians, averages and modes will differ. For most relevant 
parameters in this study the differences between the median 
and the average are small. We do not consider using a mode, 
because due to the discreteness of the model parameter space, 
the mode can be stochastically offset from the bulk of the 
PDF, and is sensitive to PDF binning. For similar reasons 
we do not use a single model with the best x^ (i-S-- the maxi- 
mum likelihood model) as representative of parameter values, 
but we do keep track of the best x^ values in order to evaluate 
the overall quality of the fitting for a given galaxy 04.31 1. Ad- 
ditionally, for certain parameters we preserve more detailed 
information on the shape of the PDF (in the case of the spe- 
cific SFR and stellar mass, we actually keep their mutual two 
dimensional PDFs). In certain cases we will discuss the er- 
rors of the derived parameters. We use 1/4 of the 2.5-97.5 
percentile range as a proxy for what would have been a Icr 
error in the gaussian distribution. We will refer to these quan- 
tities as "formal" errors. 

We prefer the abov e-described Bayesian approach (e.g., 
[Leonard & Hsu I l2003h to often-used maximum likelihood, 
since it allows one to transparently test the dependence on 
the prior, i.e., how well the data constrain certain parame- 
ters given that the shape of the resulting PDF can to some 
extent depend on the distribution of the given parameter in the 
model library (especially when the observational constraints 
are poor). Another advantage is that while the maximum- 
likelihood model might indeed correspond to a high proba- 
bility density, in the case of extensive models such as ours, it 
might actually correspond to a negligible probability "mass", 
and might therefore lead to misle ading parameter estimation 
(see for instance the discussion in lMacKay1l2003l) . Note also 
that our technique properly accounts for the degeneracies in- 
herent in some galaxy properties. Imagine that we have two 
galaxy parameters whose values can be picked in such a way 
as to produce identical model SEDs. In such case our fitting 
will give equal probabilities to all of the degenerate models, 
and the resulting PDFs for each of these two parameters will 
therefore reflect the entire range that these parameters can take 
(i.e., there will be nothing to constrain the PDFs). If we then 
characterize the errors of those parameters using the width of 
the PDF, we will obtain very large values, i.e., we will know 
that a certain parameter is poorly constrained. This, in effect, 
is in contrast to fitting using a single best-fit model, which 
would pick one of the degenerate models as the best-fitting 
(and the corresponding parameter value), without giving a re- 
searcher an idea that there is an underlying degeneracy mak- 
ing the obtained value to be quite uncertain. Finally, we want 
to emphasize that the widths of the PDFs (and also their typ- 
ical values) have a meaning only if the observational errors 
have been estimated reasonably correctly, and if the overall fit 



(which can be well characterized by the x^ of the best-fitting 
model) is reasonably good ( ^4.3l l. 

4.3. Quality of a fit 

For our initial sample of 49,346 galaxies we obtain model 
parameters in the way described in 34.21 For each galaxy we 
also note the x^ of the single best-fitting model, Xbest- ^^is 
value serves as an indicator of the overall quality of the fitting 
for a given galaxy. Since we include galaxies close to the 
edge of GALEX field of view where we know that image PSFs 
get distorted, we would first like to check if this affects the 
photometry, and therefore the quaUty of SED fitting. Thus 
we look at the average Xbast ^ function of distance from the 
GALEX image center. We find no correlation with the distance 
or degradation close to the edge. 

In the ideal case, Xbest values would form a theoretical x^ 
distribution. The average of such distribution should equal 
the degrees of freedom (in our case 6). Before looking at 
the x^ distribution we notice that Xbest values are systemati- 
cally higher for blue galaxies compared to the red. We believe 
that this difference arises from the presence of emission lines 
in blue galaxies that effectively increase the discrepancy be- 
tween the observed SED and the model (line-free) continuum. 
Therefore, we start by comparing the observed and the theo- 
retical x^ distribution of red galaxies. Initially, the observed 
distribution was shifted towards the values that were too low. 
This was indicative of photometry errors being overestimated. 
We arrive at the x^ distribution that was well-matched to the 
theoretical one when we adopt 0.01 mag calibration errors in 
each of the 5 SDSS bands."' While now the bulk of the red 
galaxies follows a x^-like distribution, there is also present a 
tail of high Xbest values, in numbers well above any expected 
distribution. One can imagine that the main reason for the 
presence of x^ outliers is that the data do not reflect the actual 
SED for some reason. We will investigate possible causes mo- 
mentarily. The high-x^ tail of red galaxies contains some 2% 
of objects. We decide to exclude these galaxies from further 
analysis (although we do account for them statistically when 
making any kind of space density calculation). Assuming that 
the reasons driving extreme x^ values are color independent, 
we also clip 2% of the blue-galaxy Xbest distribution tail. This 
leaves us with 48,295 galaxies — our final sample. 

While removing the outliers is necessary to have results not 
be biased by unreliable values, it is instructive to understand 
what causes some x^ to assume such high values. We thus 
examine 22 galaxies with extreme Xbest values, distributed 
evenly in the -1 < NUV-r < 7 range (thus sampling galaxies 
with various SF histories). In 10 cases we find obvious prob- 
lems with cataloged photometry. This is most often the case 
with SDSS (9 objects), where the visual inspection of broad- 
band SEDs and comparison with spectra indicates that some 
flux points are large outliers. Next, in 7 cases we are deal- 
ing with either close galaxies being blended in GALEX (4), or 
with SDSS shredding galaxies into smaller photometric ob- 
jects (2). For one object it is unclear whether we have blend- 
ing or shredding. Further, in one case an object shows a clear 
QSO spectrum, i.e., the continuum is affected by the AGN, 
yet it is spectroscopically (mis)classified as a galaxy. Finally, 

The "official" SDSS calibration errors are given as 2% for r band, and 
between 2% and 3% for each of the four colors. However, the individual 
magnitude errors (which figure in Equation [T) are correlated, which is why 
we find them to be smaller than 2%. Our results stay practically the same 
even if errors are increased, but the x~ distribution then deviates from the 
expected one. 
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in four cases it appears that the broad-band SEDs (and spec- 
tra) are genuinely unusual. While their number is too small to 
affect the present study, this type of outliers deserves a spe- 
cial scrutiny. It needs to be explored if there exists any set of 
model input parameters (SFH, dust attenuation, metallicity) 
beyond our already wide range, which can produce a model 
SED matching these observations. For our full sample we also 
expect on the order of 10 cases where a foreground star of sim- 
ilar brightness is superimposed with a galaxy, and falls within 
the fiber, thus producing a "composite" galaxy/star spectrum 
which would fail the fitting. 

Since our model libraries are calculated at some fixed red- 
shifts, there is some error associated with the library redshift 
not perfectly matching the galaxy redshift. We thus contrast 
the Xbest values of galaxies having the maximum redshift off- 
set (close to 0.025) to those with no offset. We find a system- 
atic increase in xlest °f 10%, which is too low to affect the 
overall distribution. 

Finally, we want to assess the effects of the "resolution" of 
our model parameters, i.e, we can ask how extensive should a 
model library be that the resolution is not an issue. To that ex- 
tent, we artificially degrade our model libraries by running full 
SED fitting using every 2nd, every 4th and every 8th model 
SED. We then compare average Xbest values from each run. 
Overall, the quality of the fits are surprisingly stable. The run 
with every 2nd model produces only a 4% increase in Xb^^^ 
values, while the run using every 16th model is 16% worse 
compared to the full library. This implies that increasing the 
extent of the library (for the given range of input parameters) 
would not bring significant improvements. This, of course, is 
not a general conclusion, but rather states that with the pre- 
cision of our data the current number of models is adequate, 
i.e., having a finer resolution of input model parameters would 
create model SEDs that are degenerate from the observational 
point of view. In addition to this, the Bruzu al & Charloti 
^2003) models themselves are not likely to be of such an ac- 
curacy to warrant more extensive model library. 

4.4. K corrections and volume corrections 

For deriving the galaxy properties from the SED fitting, we 
are not required to know the K corrections. This is because 
for a galaxy at a given redshift we compare the observed 
magnitudes to model magnitudes with bandpasses that have 
been shifted to the matching galaxy's redshift (i.e., the red- 
shift of the model library). Therefore, the "K corrections" 
are present in the SED fitting implicitly. However, for some 
applications, such a s co nstructing a CMD 03.3I ). or obtain- 
ing Vmax weights (§ 37l8b . we need to know the K corrections 
explicitly. The process of SED modeling allows us in prin- 
ciple to derive the K corrections alongside galaxy parame- 
ters. It even has the advantage of allowing one to construct 
a PDF for any K-corrected magnitude, which thus yields an 
estimate of a K correction error. We find this approach to 
be especially useful in K correcting higher redshift (z ~ 1) 
samples (Salim et al., in prep.) In the current case, given 
that our redshift range is covered by only 5 model libraries, 
we choose to calculate the K corrections using the pubhcly 
available code KCORRECT V4_l_4'^ (iBlanton & RoweisI 
I2007h . The code allows GALEX magnitudes to be used to- 
gether with SDSS magnitudes to constrain the SED fit from 
which K corrections are derived. We obtain K corrections us- 
ing the same combined GALEX/SDSS photometry used for 

" Available from http : / / cosmo . nyu . eciu/blanton/kcorrect| 



the SED fitting, except that < 3(7 UV detections are not used. 
We derive K corrections without a priori assuming an evolu- 
tion of the luminosity function. We will explicitly correct for 
evolution where appropriate. 

In cases when we require a volume-corrected sample, we 
assign a weight to each galaxy according to its Vmax value — 
the volume in which a galaxy would be visible taking into 
account redshift and apparent magnitude limits, and the solid 
angle of a survey. While our faint limit (r = 17.77) is also 
the nominal limit of the SDSS main spectroscopic survey, we 
take into account that the actual spectroscopic magnitude limit 
varies from one spectroscopic plate to another (i.e., that it can 
be brighter than r = 17.77). Thus, for each object we mod- 
ify the volume by the spectroscopic completeness (usually 
around 0.9). We take spectroscopic lim its and co mpleteness 
values from NYU-VAGC catalog (Bla nton et al.ll2005 ). We 
caution that in order to calculate Vmax correctly, it has to be 
done iteratively, because the K correction depends on redshift 
limits that determine Vmax- 

5. COMPARISON OF "UV" AND "Ha" STAR FORMATION RATES 

5.1. "UV" SFR estimates 

We derive our SFRs from photometric constraints that ex- 
tend from the UV to the z-band. However, it is the UV lu- 
minosity that constrains the SFR the most. Therefore, in the 
remainder of the paper we will refer to SFRs derived from 
the SED fitting as "UV", where quotation marks stand to re- 
mind us that it is not only UV information that contributed to 
these SFR estimates. In our model SED libraries we report 
"current" star formation rates averaged over several time in- 
tervals, most notably over the most recent 10 Myr, 100 Myr 
and 1 Gyr. By comparing the formal errors of SFR estimates 
(derived from 95 percentile range of each galaxy's SFR PDF), 
we find that the best-constrained SFR for the overall sample 
is the one averaged over the last 100 Myr (Table [TJ, which 
is the timescale for UV bright stars. For the star-forming 
class specifically, we find that the SFRs have slightly smaller 
formal errors over 1 Gyr timescales than over 100 Myr On 
the other hand, SFR estimates on timescale of 10 Myr for SF 
galaxies (a timescale that matches Ha SFR) have drastically 
larger formal errors, as expected since we do not have infor- 
mation to constrain them. In the remainder of the paper we 
will use SFR estimates averaged over 100 Myr, which repre- 
sent a good compromise between the quality and a timescale 
that is not too long compared to that of the Ha SFRs. 

5.2. "Ha" SFR estimates 

B04 have developed a method of deriving SFRs from SDSS 
spectra that is primarily based on the intensity of an Ha 
line. This represents the largest currently available sam- 
ple of precise Ha-based SFRs. Obtaining accurate SFRs 
from SDSS spectra is far from straightforward. For the 
purposes of this paper we briefly describe the procedures 
employed in B04. The essence of the B04 approach was 
to model SDSS spectra by first removing the absorption 
line spectrum usin g a combination of burst models from 
iBruzual & Chari ot (2003) , and then to model the emission 
lines using the Chariot & L onghettil (12001.) models which 
combine the Bruzual & Charloi (ll993h galaxy evolution mod - 
els with the CLOUDY photoionization code (lFerlandl ll996'). 
Al so, the same two-component dust attenuation prescription 
of ICharlot & l?all (2000) was used by B04 to model attenua- 
tion of the emission lines, and by us to model broad-band con- 
tinuum attenuation. While conventionally the attenuation of 
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Ha flux is determined by comparing the observed Ha to H/3 
ratio (Balmer decrement) to the theoretical one, the attenua- 
tion estimate in B04 is constrained using many emission lines 
(although it is dominated by Ha/H/3). Using a suite of mod- 
els, B04 apply a Bayesian approach to produce the probability 
distributions for each of the four parameters (gas metallicity, 
ionization parameter, dust attenuation, and dust-to-metal ra- 
tio). B04 thus simultaneously produce attenuation estimate as 
well as the attenuation-corrected SFR within SDSS 3" fiber. 
SFR within the fiber is constrained by many emission lines, 
with the greatest weight carried by Ha. This is one reason 
why we will denote B04 SFR and attenuation estimates by 
"Ha", with quotation marks again serving to indicate that not 
only Ha was involved in those estimates. The above proce- 
dure is directly applicable only to galaxies for which the ion- 
izing source is predominantly star formation, and not, for ex- 
ample, an AGN activity. B04 therefore calculate fiber SFRs 
directly from emission lines only for galaxies classified in the 
BPT diagram as star forming and low-S/N star forming 03. 2b . 
For other classes of galaxies (AGN, composite, galaxies with 
no Ha) for which either there is no emission line detection or 
the lines are contaminated by a non-SF ionizing source, B04 
use the relation between the D4000 spectral index and the spe- 
cific SFR (SFR normalized by stellar mass) that has been cal- 
ibrated using the star-forming galaxies. Therefore, for these 
galaxies, emission-line SFRs are used indirectly. This is an- 
other reason why we denote B04 SFRs as "Ha" w ith q uota- 
tion marks. We will return to this calibration in ^5.31 The 
above procedure gives only the SFRs within the 3" aperture 
of SDSS spectroscopic fibers, i.e, the fiber SFRs. In order to 
obtain total SFRs, an aperture correction is required. The total 
SFR is the sum of the SFR within the fiber and outside of it. 
To estimate SFR outside of the fiber, galaxies are first divided 
into a color-color grid, based on their fiber colors. Then, for 
each color cell, a distribution of / band-luminosity normalized 
SFRs (SFR/L,) is constructed by adding up all galaxies of SF 
or low-S/N SF class that fall within the given color cell. Then, 
the color of a galaxy outside of fiber is used to select the ap- 
propriate color cell, and the SFR/L, in that cell is used with 
the /-band luminosity to derive the SFR outside of fiber This 
is added to previously found SFR within the fiber to obtain the 
total SFR. In other words, the aperture correction procedure 
is based on two assumptions. First, that SFRs within fibers 
have the same dependency on color as the SFR outside of it, 
and second, that this dependency, calibrated with star-forming 
galaxies, applies to other galaxy classes. The aperture correc- 
tion factors range from close to one (no correction) to around 
a hundred-fold. On average they are 0.9 dex for the entire 
sample and 0.6 dex for the star-forming class. 

5.3. Comparison of "UV" and "Ha" SFRs for all galaxies 

The comparison between SF indicators serves the obvious 
purpose of providing a better understanding of each indica- 
tor, but also ensures the mutual check on the reliability of the 
techniques used to produce the SFR estimates — in this case, 
the SFRs based on UV (this study) and on Ha (from B04). 

We begin by comparing the formal errors of the two meth- 
ods. Note that B04 used a sample based on the expanded 
SDSS DRl. Here, when we refer to B04 SFRs we mean the 
B04 calculations applied to SDSS DR4, and available as a part 
of MPA/JHU SDSS value-added catalog. In Table[I]we show 
average formal errors in SFR for the entire sample, and for 
each galaxy class. We see that the formal errors of the two SF 
estimates happen to be comparable. 





2 


All / 




1 








•X. 















log SFR( 


-1 










-2 














-3 


'/< 1 1 1 1 " 



-3-2-10 1 2 
log SFR("UV") 

Fig. 2. — Compai'ison of dust-corrected SFRs derived bv lBrinchmann et^ 
(2004) (B04, "Ha") and the dust-corrected SFRs from this study ("UV"). 
Comparison is given for all galaxies in the sample. Our ("UV") SFRs come 
from the modeling of the broad-band SED (two GALEX UV bands and five 
SDSS optical bands), and span almost 5 orders of magnitude. While most 
galaxies compare reasonably well, those with low "UV" SFRs are quite offset 
from the equality fine. The error bar represents average formal err ors (from 
Bayesian fitting) of the two estimates. Both estimates are given for jChabrieil 
.2Q0a) IMF 
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Fig. 3. — Comparison of B04 ("Ha") and "UV" dust-corrected SFRs for 
low-S/N star-forming, SF/AGN composite, AGN, and galaxies without Ha 
detection. Note that except for low-S/N SF class, the B04 estimate comes 
from relations calibrated using the "Ha" SFRs of galaxies classified as star 
forming. In any class the discrepancies are particularly large when the "UV" 
estimate of SFR is low. Error bars represents average errors in each class. 

To faciUtate the c omparison, we f irst conv ert B04 SFRs that 
were calculated for 'Kroupal (1200 lb IMF to Chabrieii (|2003|) 
IMF used in this st i idy (b oth with 0.1-100 limits). From 
iBruzual & Charloi (l2003h models we find the conversion fac- 
tor of 1 .06 (in the sense that Kroupa IMF SFRs are slightly 
higher). In Figure|2]we plot B04 "Ha" SFR estimates against 
the "UV" for all galaxies in the sample (throughout the pa- 
per SFR always means dust-corrected SFR). Note that the 
"UV" SFRs exhibit 5 orders of magnitude of a dynamic 
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range — from all but negligible rates to those approaching 100 
Moyr"'. We notice that a large number of galaxies (those 
with higher "UV" SFRs) has comparable SFRs. However, 
there is a plume of galaxies for which B04 SFRs are up to 
2 orders of magnitude higher than the "UV" ones. This sys- 
tematic discrepancy cannot be explained by random errors of 
the two methods. Recall that B04 do not derive SFRs di- 
rectly for classes other than star-forming (and low S/N star- 
forming), but instead use relations calibrated on SF class sam- 
ple. Thus, in Figure [3] we break the comparison into galaxy 
classes (omitting for the moment the SF class). For the low- 
S/N star-forming class (upper left panel), for which B04 SFR 
estimates do come from the emission lines, the comparison 
is generally good for the majority of galaxies (apart from the 
overall offset that makes "Ha" SFRs higher). For the other 
three classes in Figure[3j B04 SFRs rely on calibrations made 
using the SF galaxies. Composite (SF/AGN) galaxies com- 
pare relatively well, but the comparison becomes significantly 
worse for AGN, where we can see both the significant offset 
for galaxies with higher SFRs, and the prominent plume for 
those with lower "UV" SFRs. Finally, for galaxies for which 
no Ha detections were possible, the discrepancy is very large 
for almost all galaxies. It is this class that contributes the most 
to the plume seen in the comparison of all galaxies (Figure 
121). Given that the galaxies without Ha detections belong al- 
most exclusively in the red sequence (Figure[T]i, the low SFRs 
rates estimat ed by the SED fitting ("UV") appear more realis- 
tic (see also 37.1b . Also, such high SFRs as estimated by B04 
would result in a high fraction of No Ha galaxies to have a 
detectable Ha emission, which is obviously not the case. In 
general, it appears that the largest discrepancies between the 
"Ha" and "UV" rates occur for galaxies with low SFRs, es- 
pecially when B04 estimate SFRs based on calibrations that 
employed SF galaxies. 

B04 SFRs are the sum of SFRs within and outside of fiber. 
To trace the source of the SFR discrepancies, we revisit the 
calibrations used to derive B04 SFRs for non-SF classes. 
First, within the fiber, B04 base the SFRs for non-SF classes 
on the relation between the specific SFR (in the fiber) and the 
D4000 index calibrated using the SF and low-S/N SF galax- 
ies. We show this calibration in both panels of Figure|4]as the 
dashed line (this can be compared to their Figure 11). Then, 
B04 use D4000 to determine the specific SFR (and therefore 
the SFR itself) of galaxies of other classes. In the upper panel 
of FigurelH we compare B04 relationship to our ("UV") rela- 
tionships for the star forming (bold solid line) and the low-S/N 
star-forming (thin solid line) galaxies. Note that in principle 
we cannot compare the B04 specific SFR inside of fiber to 
the total specific SFRs obtained from the UV, so the compar- 
ison is more for illustrative purposes. At D4000 < 1.7 either 
SF or low-S/N SF lines from "UV" compares well to B04. 
However, notice that there are no galaxies classified as star 
forming beyond D4000 = 1.9. Therefore, from that point on, 
B04 relationship is based only on low-S/N SF class which 
are most likely just AGN contaminants (see 33.2| i. In con- 
trast to B04, the "UV" relationship for low S/N SF continues 
to decline as D4000 increases. In the lower panel of Figure 
ID we show the same B04 calibration against "UV" rates for 
the AGN (thin solid line) and galaxies with no Ha detections 
(bold solid line). In these cases, the relationship for these two 
types as indicated by the "UV" is entirely different from that 
of the star forming galaxies, even at lower D4000 values. For 
clarity we omit showing the composite class. Altogether, this 
leads to a conclusion that the relationship between the specific 
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Fig. 4.— Relationships between SFR/M, (the specific SFR) and the D4000 
spectral index. In the upper panel we compare the relationship from B04 
(constructed from SF and low-S/N SF galaxies combined ) to our "UV" rela- 
tionships for the SF and the low-S/N SF class, separately. The B()4 relation- 
ship comes from measurements within the fiber, while others are total, so this 
comparison is primary an illustration. B04 use such relationship to determine 
fiber SFRs of non-SF classes of galaxies. Note that beyond D4000 = 1.9, the 
B04 calibration is actually based on low S/N spectra which may be contami- 
nated by AGN. In the lower panel the same B04 relationship is compared to 
our "UV" rates for the AGN, and for galaxies with no Hq detection. The two 
follow different relationships, with AGN and "No Ha" classes having much 
lower specific SFRs than the SF class. The non-unique mapping between the 
specific SFR and the D4000 index for various galaxy groups, and the fact that 
the congelation cannot be established well using the emission lines indicates 
that fiber SFRs in B04 may be systematically affected. 
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Fig. 5. — Relationship between the luminosity-normalized SFR and op- 
tical colors. B04 use luminosity-normalized SFRs binned in g-r vs. r-i 
color-color grid to estimate SFR outside of fiber. They build their relation- 
ships using SF galaxies, but apply them to other types as well. Here, based on 
"UV" measurements, we show that different classes of galaxies follow dif- 
ferent relationships, and that applying the one from SF galaxies would lead 
to systematic differences, which cause most of the SFR discrepancies seen in 
Figure[3] 
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SFR and the D4000 index is not unique for different classes of 
galaxies. We tested if this is an artifact of the fact that D4000 
comes from the central 3", while the specific SFR for "UV" 
is integral. We find that the relationships are qualitatively the 
same for z < 0.07 and for z > 0.12 samples for which D4000 
probes different physical sizes. Instead, the most likely expla- 
nation is that since D4000 and specific SFR are sensitive to 
SF activity on different timescales, the galaxies with differing 
SF histories will have different relationships. To much lesser 
extent, some of the differe nces could be due to metallicity 
dPo^panti & Barbarolll997h . 

Given the level of aperture corrections (a factor of several), 
the analysis given above does not provide the explanation 
for the large discrepancies we see when comparing the to- 
tal SFRs. Since the SFR outside of fiber dominates in B04 
total SFR, we will now focus on the calibration B04 use to 
determine it. As outlined in 35.2I B04 aperture correction re- 
lies on the calibration of the luminosity-normalized SFR as a 
function of g-r and r-i color. Using our UV-based SFRs 
we can test one of the two assumptions behind B04 aper- 
ture corrections — whether the SFR/L calibration against color 
holds for different classes of galaxies. In Figure |5] we plot 
SFR/L against two colors that B04 use to determine SFRs 
outside of fiber. The bold line represents the relation from 
galaxies classified as star-forming. This is the basis for B04 
calibration which they apply to other classes as well. The 
dashed and the thin solid lines show the relations for AGN 
and galaxies with no Ha. We see (left panel) that as the g-r 
color increases, the discrepancy between SF and other classes 
rises, and reaches some 0.8 dex for AGN and 1.3 dex for No 
Ha. Similar levels of difference are present against r-i (right 
panel). The level of these differences matches the discrepan- 
cies (plumes) in the comparison of B04 and "UV" SFRs in 
Figure [3] We can interpret these differences as reflecting the 
fact that the colors used for B04 calibration trace old popula- 
tions, while the current SFR traces young populations, and it 
is not surprising that the two will differ for different classes 
of galaxies. Also note that the discrepancies of specific SFRs 
of different classes of galaxies are qualitatively similar in the 
case of D4000 and in the case of g-r color, which is not sur- 
prising since they are both sensitive to population age on the 
similar timescales, and are similarly not affected by the dust. 
To conclude, UV-derived SFR is subject to fewer limitations, 
so it can be applied to a more diverse types of normal galaxies. 

5.4. Comparison of "UV" and "Ha " SFRs for star-forming 

galaxies 

To avoid the problems indicated in the previous section, we 
need to compare the SFR estimates from the "Ha" and the 
"UV" methods in galaxies where Ha is well-detected (S/N 
> 3) and arises predominantly from star formation, i.e., to 
galaxies classified as star-forming (SF). In Figure|6]we show 
the B04 "Ha" star formation rates plotted against our "UV". 
We see that the comparison, spanning some three orders of 
magnitude in SFR, is exceptionally good. Formal error bars 
are comparable, with the "UV" being somewhat smaller (see 
Table [T]i- The scatter (standard deviation of the difference) 
of the two measurements is 0.50 dex. When 3(t outliers are 
excluded, the scatter is reduced to 0.36 dex. This is very well 
matched by the sum (in quadrature) of the formal errors of 
the two methods (0.35 dex), confirming that the two measures 
are predominantly independent. There is an average offset 
between the two SFRs of only 0.06 dex in the sense that "Ha" 
SFR is higher (which reduces to 0.02 dex, i.e., 5%, when 3(7 




log SFR("UV") 

Fig. 6. — Comparison of B04 ("Ho") to our "UV" SFRs (both indepen- 
dently dust-corrected) for galaxies classified as star-forming. The two com- 
pare very well on one-to-one basis. Also, the scatter is compatible with each 
measurement's errors. 
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Fig. 7. — The difference of B04 ("Ha") and our "UV" SFRs, for galaxies 
classified as star-forming, as a function of aperture correction applied to B04 
SFRs. There is no apparent trend, suggesting that the aperture corrections 
derived by B04 are quite robust. Running average (in 0.05 dex bins) is shown 
as a thick line, and the ilcr range as dashed lines. 

outliers are excluded). 

However, Figure |6] alone can potentially hide some system- 
atic trends between the two SFR estimates. Where can the 
differences in the SFR estimates arise from? Deriving SFRs 
for either methods incorporates several steps: (a) obtaining 
full "Ha" or "UV" observed luminosity (which in the case 
of "Ha" involves aperture corrections), (b) correcting the ob- 
served luminosity for dust attenuation, and (c) converting the 
luminosity into a SFR. While both methods perform these 
three steps simultaneously, we can still design tests that would 
differentiate between. Given in some cases the large aperture 
correction applied by B04, we first plot in Figure |7] the dif- 
ference in SFR estimates (SFR residuals) with respect to the 
level of aperture correction applied to "Ha" rates. We see 
no correlation. We also check, but do not find a correlation 
between the residuals and the apparent sizes of galaxies (not 
shown). Those are good indications that the aperture correc- 
tions of B04 are quite reliable for galaxies classified as star- 
forming. Finally, we find that the SFR residuals exist when 
plotted against the stellar mass of a galaxy (shown in Figure 
[8]). We have a change of SFR residual of 0.38 dex over the 
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Fig. 8. — The difference of B04 ("Ha") and our "UV" SFRs, for galaxies 
classified as star-forming, as a function of galaxy stellar mass. There is a 
clear trend with respect to mass, leading to a difference of SFRs of 0.38 dex 
(a factor of 2.4) over the 8.5 < logM» < 1 1 range. Running average (in 0.2 
dex bins) is shown as a thick line, and the itlo" range as dashed lines. 



8.5 < log M^, < 1 1 range, with 1 cr dispersion of 0.33-0.41 dex 
around the running mean. To double check if the residuals are 
in any way connected with B04 aperture corrections, we limit 
the sample to farther, and therefore on average smaller galax- 
ies (z > 0.12). We still find that the residuals correlate with 
stellar mass. 

Before ruling out step (a) as being responsible for the mass- 
dependent residuals, we perform two additional checks. First, 
imagine there is an offset in zero points between GALEX and 
SDSS magnitudes (after all, they come from different surveys 
and are not measured in identical apertures). E.g., if GALEX 
fluxes were systematically overestimated, this could cause 
"UV" SFRs to be overestimated as well, possibly in such a 
way to preferentially boost blue, low-mass galaxies (as the 
trend in Figure [8] would suggest). To test this, we perform full 
SED-fitting runs in which we make FUV and NUV fainter 
by 0.1, 0.2 and 0.4 mag. We first check the quality of these 
new fits by comparing the distribution of xle-n values with the 
nominal (no magnitude offset) run. Runs with 0.2 and 0.4 
mag offsets produce evidently inferior fits, akeady suggesting 
that any potential offset cannot be that large. The run with 0. 1 
mag offset, however, appears as good as the nominal run, with 
an even slightly smaller average Xbest- However, comparing 
the "UV" SFRs produced with offset GALEX magnitudes to 
"Ha" SFRs, we find that while slightly flattening the slope of 
the residuals with respect to mass, the trend is not eliminated. 
In fact, the scatter of the residuals around the mean is larger 
when modified "UV" SFRs are used instead of nominal. We 
also look at the quality of fits with GALEX magnitudes offset 
in the opposite direction (making them brighter), but such fits 
are evidently inferior. 

The second test concerns SDSS magnitudes. Namely, B04 
use CMODEL SDSS magnitudes to transform fiber SFRs 
into total SFRs, while we use MODEL SDSS magnitudes 
when performing SED fitting.'^ Therefore, we perform an- 
other SED fitting using CMODEL magnitudes instead. First, 

SDSS MODEL magnitudes (defined as either the exponential or 
de Vaucouleurs magnitude, depending on which profile better describes a 
galaxy) are preferred magnitudes for SED fitting as they preserve relative 
fluxes (colors) better. CMODEL magnitudes (defined as a composite of ex- 
ponential and de Vaucouleurs magnitude) should provide a good measure of 
a total galaxy light. 
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Fig. 9. — Comparison of B04 ("Ha") to our ("UV") estimates of V-band 
dust opacity (attenuation), for the SF galaxy class. The formal error of both 
estimates (average errors shown) is large, leading to a large scatter. However, 
there is an agreement between the two in the general sense. B04 tv is to first 
order constrained by the Ha/H/3 ratio (Balmer decrement), while our estimate 
is predominately constrained by the UV slope. Both estimates were made in 
accordance with the .Chariot & Fall l ,200Q) two-component dust attenuation 
model. 
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Fig. 10.— The difference of B04 ("Ha") and our "UV" estimates of V- 
band dust opacity as a function of galaxy stellar mass (for SF galaxies). There 
is a trend with respect to mass, leading to a 0.5 difference over the 8.5 < 
log M, < 11 range. Running average (in 0.2 dex bins) is shown as a thick 
fine, and the itlo" range as dashed lines. 

we notice that the quality of fitting with CMODEL magni- 
tudes is noticeably inferior, stressing their inadequacy for pro- 
ducing reliable color estimates. Also, the use of CMODEL 
magnitudes does not remove the trend of SFR residuals, and 
the scatter of the residuals becomes larger than in the nominal 
run, especially at higher masses. With this test we exhaust the 
possibilities that discrepancies arise in step (a) above. 

We now move onto step (b), i.e., correcting the observed 
flux for dus t attenuation. Note that both B04 and this study 
use ICharlo t & Fall (2000) prescription for dust attenuation, 
the aim of which was to produce consistent treatment for Ha, 
and UV continuum attenuation. Is the presence of the SFR 
residuals an indication that this dust attenuation model does 
not produce fully consistent answers? In our SED fitting we 
keep track of estimates on Ty — the dust opacity in rest-frame 
V band. We can thus compare our ry values with those ob- 
tained by B04 (which we denote tv("UV") and ri/("Ha"), re- 
spectively). While we use the same model to constrain at- 
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Fig. 1 1 .— The difference of B04 ("Ha") and the modified "UV" SFRs as a 
function of galaxy stellar mass. We modify one of the two SFRs (in this case 
the "UV" SFR) by the difference in B04 and our attenuation estimate. After 
this moditication, the trend of the residuals with respect to mass is ehminated, 
and the overall difference is close to zero. Therefore, we find that emission 
lines and the continuum essentially produce the same SFR estimates. Run- 
ning average (in 0.2 dex bins) is shown as a thick line, and the ± lo" range as 
dashed lines. 

tenuations, B04 obtain them from the emission lines (to first 
order from the Balmer decrement), while we obtain them 
from the broad-band SED (to first order from the UV spec- 
tral slope, (3, or equivalently, the UV color). Also, B04 mea- 
surements are restricted to the fraction of the galaxy inside 
of the fiber, which can in principle produce some systematic 
differences. Direct one-to-one comparison, shown in Figure 
|9l shows a rough agreement, albeit with a large scatter due 
to relatively large errors of both estimates (average formal er- 
rors shown as the error bar). Looking instead at the distribu- 
tions of the two Tv estimates (not shown), we find them to 
be very similar on the whole, with the Ty("UV") slightly off- 
set towards the larger values with respect to TvC'Ha"). Since 
we find SFR residuals to correlate with the mass, it is more 
instructive to check whether such a trend exists in the dif- 
ference of Ty estimates from "Ha" and "UV". In Figure [TO] 
we plot Tv residuals against the stellar mass. Despite the 
large scatter, we find a systematic trend. The gradient is 0.47 
(equivalent to 0.52 mag) over the 8.5 < logM, < 11 range. 
For lower masses, "UV" Ty attenuation is higher than "Ha", 
which implies that the observed UV flux is corrected more 
than Ty("Ha") would require, resulting in "UV" SFRs be- 
ing higher than SFR("Ha"). The situation is reversed at high 
masses. This trend matches the sense of the trend of SFR 
residuals. 

In order to establish whether the difference in ry estimates 
is the dominant cause of the trend of SFR residuals, we will 
try to modify one of the two SFR estimates (e.g., the "UV" 
SFR) based on the difference in Ty values. Note that this 
is an approximate technique, since our SED SFRs are not 
"corrected" for dust by the applicatio n of a single number , 
but rather by a complex application of ICharlot & Falll ( l2000h 
prescription on many different populations that constitute the 
model SED. Nevertheless, for the purposes of this test we 
would assume that the attenuation in the FUV drives the cor- 
rection of the "UV" SFR estimate. We therefore modify it 
with: 

SFR("UV")mod = SFR("UV") x lo''^'^'*™'' , (3) 
where AAf (yy is the difference between the attenuation in the 



rest-frame FUV as implied by Tv("Ha") and that determined 
from the SED fitting. We look at the parameters from the SED 
fitting to calibrate the relationship between the attenuation in 
FUV and the V-band opacity, and find the following relation- 
ship: 



AAff/v = /i(5.4-0.84Tv)Ary, 



(4) 



where Ary = ry("Ha") - tv("UV"), while ry and /i come 
from the SED fitting (i.e., they are "U V" measurements). The 
parameter ji is the coefficient in the ICharlot & Falll (l2000l) 
model that determines the fraction of the total opacity af- 
fecting the diffuse ISM. We modify SFR("UV") according to 
Equations [3] and m and in Figure [TT] re-plot the difference in 
SFRs against the mass. The trend of the SFR residuals with 
mass has disappeared, and the difference stays within ±10% 
across the full mass range. The dispersion around the running 
mean is 20% larger than with the non-modified SFR resid- 
uals, which is to be expected due to the approximate nature 
of the modification applied. Therefore, we conclude that the 
"UV" continuum and the "Ha" in essence produce identical 
answers for the star formation rate, the apparent difference 
stemming from different estimates of the attenuation that the 
two methods provide. With this in mind, we can also rule out 
other possibilities for the cause of the original differences in 
the SFRs, such as the reliability of converting the luminosity 
into a SFR (point (c) abo ve). Also, these results sh ow that, 
contrary to some notions, iBruzual & Chariot (l2003h models 
cannot be too much off in the UV. 

Although we can now account for the SFR differences, we 
would still like to understand the origin of attenuation differ- 
ences that cause them. Related to this, we should try to de- 
termine which attenuation estimate (emission-line or the UV 
continuum) is more accurate. At this point we need to em- 
phasize that our SED fitting has a non-flat prior distribution 
of Ty (described in 34.11 ). while B04 uses a flat prior, with at- 
tenuation taking values in the 0.01 < Ty < 4 range. While the 
non-flat prior used in our study is a reasonable assumption, 
and the one that helps constrain the parameter space to physi- 
cally realistic values, we would like to check whether the dif- 
ferences in priors induce any systematic effects. To that effect 
we first perform an additional SED fitting run, with the model 
libraries created with a flat attenuation prior, taking values in 
the < Ty < 3 range. At lower (ry < 0.7) values, we find the 
distribution of the resulting Ty values to be similar to the origi- 
nal distribution, while at higher values it displaces the original 
peak, located at ry « 1 .3, to a broader distribution peaking at 
Ty « 2.0. The latter is the obvious consequence of assign- 
ing equal probability to models with high attenuations. In any 
case, we end up with a ry distribution that differs more with 
respect to the B04 Ty("Ha") than did the original ry("UV"). 
Nevertheless, we proceed and check the relationship between 
the new SFR residuals and the stellar mass. The trend ob- 
served in the original relationship is still present, although it is 
somewhat weaker at low masses. Anyhow, it appears that we 
cannot force the two attenuation estimates to reach an exact 
agreement by a simple modification of the prior distribution 
used in SED fitting. 

Next we will try to get a sense of which ry estimate is more 
realistic. Since the direct comparison with an external in- 
dependent measurement is not readily available, we will try 
to evaluate the differences using internal relations. For star- 
forming galaxies we expect attenuation to be correlated with 
the stellar mass. In Figure [T2l we show this relation for our 
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Fig. 12. — Attenuation-mass relationship for star-forming galaxies. The 
two attenuation estimates, one from our SED fitting ("UV", upper panel), 
and the other from B04 ("Ha", lower panel) are plotted against the stellar 
mass. The B04 estimate (principally from Ha/H/3 ratio) produces a tighter 
relationship with a greater span of attenuation values, indicating that it is 
probably more accurate. 

"UV" estimates (upper panel) and B04 "Ha" estimates (lower 
panel). It is apparent that the latter defines a more tight rela- 
tionship. Also, rv("Ha") spans a larger range of values than 
rv("UV") (1.4 vs. 0.8 over 8.5 < logM, < 11). Finally,, the 
dispersion around a running mean is some 20% smaller for 
TvCHa"). This all indicates that B04 ry estimates are prob- 
ably more accurate. While the two tv estimates agree well 
on the whole (Figure |9|l, the "UV" estimate appears more 
poorly constrained, which may cause it to more often take val- 
ues away from the extrema, suppressing the span in relation- 
ship with the mass. As a result, the attenuation estimates for 
low-mass galaxies will be overestimated, and for high-mass 
galaxies underestimated. It is difficult to directly test this ex- 
planation. Namely, if we restrict the analysis only to galaxies 
with small formal errors in tv("UV") (or alternatively small 
errors in FUV -NUV), we at the same time bias the sam- 
ple to s maller values of ti/("UV"). Alternative explanation is 
that the ICharlot & Fa ll (2000) oc A"°^ extinction curve devi- 
ates from the true extinction curve in a way that would be 
mass-dependent (or some quantity related to mass, such as the 
metallicity). A possibility that the extinction curve is mass- 
dependent is raised in Johnson et al. (2007). We would be 
able to acquire further insight by comparing both attenua- 
tion estima tes to some ex t ernal measure, such as the IR ex- 
cess (IRX). Hohnson et al.l (|2007|) have recently obtained IRX 
attenuations from Spitzer MIPS and GALEX observations of 
SDSS galaxies, therefore such sample can be used to include 
UV (i.e., SED) derived attenuations. Also, the SED-based at- 
tenuation estimates can possibly be improved by adding the 
near -IR data, such as the J HKs photometry from 2MASS sur- 
vey (Skrut skie et al.ll200^ . which would improve attenuation 
estimates by placing stronger constraints on the stellar metal- 
licity. 
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Fig. 13. — Relationship between the attenuation in the GALEX far-UV 
and the UV slope. The FUV attenuation estimate comes from our SED fit- 
ting, and the FUV —NUV rest-frame color serves as an indicator of the UV 
SED slope. Data and the solid line linear fit is for galaxies classified as star 
forming. The dashed line is a fit to blue (rest frame NUV - r < 4) galaxies. 
The fits reproduce the trends in general, however, there is a significant scat- 
ter. Meurer et al. 1 1999) relation for starburst galaxies is plotted as the dotted 
line (M99). It is obviously not appropriate for our sample of mostly normal 
star-forming galaxies. 

6. OBTAINING THE UV SFRS WITHOUT SED MODELING 

In this work we are using sophisticated SED modeling to 
obtain UV-based star formation rates. In many applications 
such a detailed approach is not practical, or not even possi- 
ble. In such cases one would like to obtain a reasonably good 
estimate of a dust-corrected star formation by applying some 
simple transformations to the UV photometry. While such 
methods have been used in many previous studies, here we 
will use the results of the detailed SED analysis to calibrate 
such simplified models, with a special emphasis on users of 
UV data obtained with GALEX. Since for the foreseeable fu- 
ture GALEX will remain the only facility capable of observing 
a large number of galaxies in the UV, it is not without justifi- 
cation to treat its filters as defining some standard photomet- 
ric bands in the UV domain. Researchers who study high-z 
galaxies in rest-frame UV can calibrate their blueshifted fil- 
ters against GALEX FW and NUV response curves.'^ 

Obtaining a UV-based SFR consists of K-correcting the UV 
magnitudes, estimating the dust attenuation of the FUV flux, 
and converting the dust-corrected FUV luminosity into a SFR. 
We assume that the user will correct their data for the Galactic 
extinction, and then apply some standard K-correction pro- 
cedure to obtain rest-frame FUV and NUV magnitudes. Of 
course, obtaining reliable K corrections requires optical or 
near-IR photometry. To calibrate the relations in this section, 
we use KCORRECT V4_l_4 ( ^4^41 1. 

We start by using the results of the full SED fitting to cali- 
brate the well-known correlation between the attenuation and 
the UV spectral slope (e.g., Calzetti et al. 1994). This is of- 
ten refereed to as IRX-/3 relation, although strictly speaking 
IRX indicates IR excess, which is correlated with UV attenu- 
ation. IRX-/9 relationship for normal galaxies, such as those 
in our sample, has been previously studied with GALE X data 
by Seibert et al. (2005); Cortese et al. (2006); Gil de Pa z et al.l 
(2007); iBoissier et all (120071) ; IPanuzzo et all (120071) . In Fig- 
ure [13] we show FUV attenuations (Afuv) obtained from our 

GALilZ filter response curves can be obtained from 

|http : //galexgi . gsf c ■ nasa . gov/tools/Resolution_Response| 
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SED modeling, for galaxies classified as star-forming, plot- 
ted against the rest-frame UV color, which is linearly corre- 
lated to the UV spectral slope. We see that the majority of 
objects lies along the ridge. This confirms that there exists 
an "IRX-/3" relationship for normal galaxies, not just for star- 
burst galaxies as usually assumed. To quantify the relation- 
ship for our sample, and thus allow Afuv to be determined 
from GALEX observations of normal galaxies, we fit a linear 
function to running medians (with value from each 0.05 mag 
bin weighted equally). Using the medians is necessary in or- 
der to avoid the fit to be affected by numerous outliers. Also, 
we find that past some red UV color the attenuation does not 
seem to rise, so we adopt a constant value. In any case, Afuv 
exhibits a large scatter for such red colors. 

^( 3.32^\FUV-NUV) + 0.22, ^^(FUV -NUV) < 0.95 
™^ ~ \ 3.37, ^^(FUV-NUV) > 0.95, 

(5) 

where Afuv is in magnitudes, and superscript designates 
rest-frame colors. This relation is plotted in Figure [13] as a 
solid line. In many applications, especially at higher red- 
shift, spectroscopic classification may not be available. In 
those cases one can select star-forming galaxies based on their 
color. For blue-sequence galaxies (•^(NUV - r) < 4, and with- 
out applying any class selection) one should use a slightly 
modified relation: 



2.99^\FUV-NUV) + 0.21, ^\FUV-NUV) < 0.90 
2.96, ^\FUV-NUV)> 0.90, 

(6) 

shown in Figure [13] as a dashed line. Note that these rela- 
tionships are optimized for GALEX bandpasses and for nor- 
mal star-forming galaxies. We over plot Meurer et al. (1999) 
relation as a dotted line (we used fSeibert et alj|2005l transfor- 
mation to obtain FUV -NUV from the UV spectral slope). 
Apparently, this relation, constructed from a sample of star- 
burst galaxies, is well above the majority of galaxies in our 
sample. Seibertet al. (2005|) and lCortese et al. (2006), using 
smaller sa mples (RC3 and clu ster galaxies, respectively) also 
find that Meu rer et al.l ( 1999h relation overpredicts FUV at- 
tenuation. Seib ert et al. I (120051) derive a relati on that exactly 
bisects our SF fit and the 'Me urer et alJ (1 19991) relation. Also 
note that these relations pertain to the optically selected sam- 
ple (since practically all SF galaxies are detected in UV bands, 
there is no additional UV selection). In general, samples se- 
lected at different wavelen gths will have diff erent IRX-/3 re- 
lations, as demonstrated by iBuat et al.l (l2005b . using GALEX 
NUV and IRAS 60/im selected samples. 

For a source with known redshift, the dust-corrected, rest- 
frame FUV is converted into a FUV luminosity. The final 
step requires converting the luminosity into SFR. Following 
B04 notation, we define: 



7fFuv = L%y/SFRCm"), 



(7) 



to be the inverse conversion factor between a dust-corrected 
rest-frame FUV luminosity (in erg s"' Hz"') and the SED- 
derived SFR. The conversion factors comes from the stellar 
population modeling that was used to perform the SED fit- 
ting. Like the equivalent conversion factor for Ha luminosity, 
Vfuv sensitive to metallicity, albeit more weakly so. For 
our sample (the metallicity of which is on average O.SZ©), 
the median conversion factor is: 
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Fig. 14. — Comparison of the SFRs "predicted" using a set of simple trans- 
formations, to SFRs derived using the full SED fitting. The predicted SFR 
is obtained using the transformations that we calibrate from the SED SFRs. 
They are applied to K-coirected FUV and NUV magnitudes (see The 
comparison is good for the majority of galaxies (greyscale scatter plot). The 
dashed lines show 16 and 84 percentile values of the conditional distribution 
at each SFR("UV"), from which we see that the predicted SFRs deviate from 
SFR("UV") at low and high SFRs. 



Iog7;0f,^ = 28.165. (8) 

The above factor is given for lChabrieii (l2003h IMF, which 
is used in this paper The often-used conversion factor given 
by Kennicutt (1998) for the Salpetej ( Il955h IMF is: 

SFR(M0yr-') = 1 .4 x lO'^^L^ieTg-h-^Uz-^). (9) 



We find using i Bruzual & Chariot ( l2003h models that the ap- 
propriate transformation factor between UV-derived SFRs 
that assume Chabrier and Salpeter IMFs is 1.58 (both with 
0.1-100 M© limits)^". Thus our "empirical" conversion factor 
(Equation[S]l for Salpeter IMF becomes: 



SFR(Moyr-')= 1.08 x 10 



-^'L%y(crg 



's-'Hz"'). (10) 

This implies that the iKennicutl (Il998h conversion factor 
(Equation O is 30% hig her. We verify th at the effect of the 
difference in bandpasses ( Kennicutt! 1 9981 factor was given for 
the 1500-2800 A range, while the GALEX FUV filter spans 
1300-1800 A) is completely negligible. The actual reasons 
for the difference are two-fold. First, even wh en obtaining the 
model conversion factor using iBruzual & Qia rlot (2003) and 
the same assumptions as given by'K ennicuttI fl998) (i.e., solar 
metalhcity and a constant star formation history), we still find 
iKennicutig (11998) conversion to be 15% hi gher (i.e., the corre- 
spond ing rj is lower than one produced by IBruzual & Chariot! 
I2003h . Further, we have that our sample has average metallic- 
ity somewhat lower than the solar. This accounts for another 
5% difference. Finally, the remaining 10% difference stems 
from the fact that our sample has a variety of star formation 
histories, giving rj that is on average higher than that for con- 
stant SF history. To summarize, for optically selected sam- 
ples similar to ours, we suggest using "empirical" conversion 
given in Equation[8]or[Tol 

Finally, we combine the above ingredients and "predict" the 
SFR from K-corrected FUV and NUV magnitudes. In Figure 

Note that the conversion between the IMFs is not generally the same for 
Ho-derived SFR, UV SFR, or stellar mass. 
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[14] we compare the SFR obtained using this simple prescrip- 
tion to the original SFR("UV")- We observe a good agreement 
for the majority of objects (the greyscale scatter plot). In addi- 
tion, for each bin in SFR("UV") we show a 16-84 percentile 
range of the predicted SFR values, i.e., their conditional distri- 
bution. This allows us to track the trend in the regions where 
there are few points to show on the greyscale scatter plot. We 
see that the predicted SFR starts to deviate at both the low 
and the high SFR ends, by up to 1 dex. Good agreement is 
obtained in the -1 < log SFR < 1 range, i.e., where the ma- 
jority of galaxies lie. The reason for the deviations at low and 
high-SFR end stem from the inability to estimate Afuv reli- 
ably only from FUV and NUV for galaxies having extreme 
SF. No simple modifications of Equations |5] and |6] are able 
to reproduce Afuv estimate of galaxies with extreme SFRs 
without destroying the agreement for the majority of normal 
galaxies. In -1 < log SFR < 1 range, the 16-84 percentile 
range is 0.35 dex wide, indicating that the predicted values 
have an equivalent la scatter of 0.17 dex. Therefore, the to- 
tal error of the predicted SFRs in this range is 0.26 dex (with 
the error of SFR("UV") itself added in quadrature). This is a 
small penalty considering the ease of the method. Note that 
the prescription given here is optimized to produce good re- 
sults for the large majority of galaxies, and that a result for 
any given galaxy should be used with caution. 

7. STAR FORMATION IN THE LOCAL UNIVERSE 

7.1. Star formation history and stellar mass 

It is now commonly accepted that the star formation history 
of a given galaxy depends strongly on its mas s (for early work 
based on UV see Boselli et al. 2001; Gav azzi & Scodeggiol 
Il996l) . and to some extent its environment. To first order, we 
can characterize the star formation history using the specific 
star formation rate, which is defined as the ratio of the current 
SFR to the current stellar mass. Thus, higher values of the 
specific SFR indicate that a larger fraction of stars was formed 
recently. The unit of specific SFR is inverse time, and it is 
thus often referred to as the galaxy build-up time. However, 
this should not be confused with an actual galaxy age, since 
the specific SFR tells us only how long it would have taken to 
build a galaxy assuming it had a current SFR throughout its 
lifetime, and assuming that the current stellar mass is close to 
the total mass formed in stars during a galaxy's lifetime (the 
two are not the same because of gas recycling). 

For each galaxy in our sample we construct a 
two-dimensional probability distribution function in 
(log(SFR/M*),logM*) plane as outlined in ^ We use 
full 2-D PDF instead of collapsing it into a single value in 
order to better represent low-SFR galaxies whose PDFs are 
often very wide (i.e., not well constrained) in specific SFR, 
and in which cases the single value, such as the average of 
the marginalized PDF, would tend to artificially concentrate 
in the middle of the PDF. In Figure[T5]we show the co-added 
2-D PDFs in 0.1/0.05 dex wide bins. In the upper panel 
of Figure [15] as in all previous figures in the paper, we 
show unweighted values, i.e., each galaxy in our sample 
contributes equally to the greyscale density. Given our 
redshift and magnitude limits, we adopt the r-band absolute 
magnitude limit of Mr = -14.75 (which we verify by con- 
structing a luminosity function). This luminosity limit leads 
to the completeness in stellar mass of logM* w 7.5, which 
is below the plotted mass range. The dashed line is shown 
for reference and represents a constant SFR of lM0yr"'. 
We see that lower-mass galaxies appear to be confined in 
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Fig. 15. — Dependence of the star formation history on the stellar mass. 
We use the specific SFR (SFR/A/, ) as an indicator of a star formation his- 
tory. Galaxies with a larger fraction of recent star formation will have a 
higher value of the specific SFR. Instead of a single value, each galaxy is 
represented with a full two dimensional probability function. The dashed 
line, shown for reference, represents a constant SFR of 1 Mq yr"' . The up- 
per panel gives equal weight to every galaxy in the sample, while the lower 
panel shows logarithms of volume-coirected values (i.e., weighted by Vmiix-) 
Uneven behavior at low masses is because of a small number of galaxies (or 
no galaxies) in some bins. Note that the sample is only optically selected. 
Our completeness hmit is below the lowest mass on the plot. 

a relatively narrow sequence which declines as the mass 
increases. At masses above 10'°Mq, some galaxies assume a 
much lower value of the specific SFR, i.e., the sequence turns 
over In the lower panel, each galaxy's PDF is weighted by 
its Vmax value, i.e., the sample is volume-corrected. Since the 
volume-corrected values span a very wide range of values, we 
now show logarithms of number densities, scaled so that the 
full greyscale range displays densities from maximum (black) 
to 2 dex below the maximum value (lightest shade of gray). 
Each horizontal slice in effect represents a mass function of 
galaxies of a given specific SFR. It now becomes apparent 
that the sequence seen in the unweighted plot extends to 
lower masses, with an increasing number density of galaxies. 
Also, the mass distributions at low specific SFR extends 
to lower masses than implied by the non-corrected plot. 
Considering the size of the sample, the accuracy of SFRs, and 
the sensitivity to low levels of SF, this figure represents the 
best available assessment of the star formation history in the 
local universe. The typical error in stellar mass is only 0.07 
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Fig. 16. — Conditional dependence o f the star formation histoiy on the 
stellar mass. Same sample as in Figure [Ts] but with each 0.05 dex wide 
mass bin normalized to its maximum separately. This allows us to see what 
specific SFRs dominate at each mass, and to follow activity where there are 
to few objects to show in the standard plot. The lower panel shows running 
modes from the upper panel together with a Schechter-like fit to those points 

dex, and 0.40 dex in the specific SFR (note that the errors in 
specific SFR vary by a factor of three across the sample). We 
complement the above figure with Figure [16] (upper panel), 
where we display conditional probability distribution in each 
0.05 dex wide mass bin, which is equivalent to normalizing 
each column to its maximum value. Such representation 
allows us to see what specific SFR dominates at a given 
mass, and it also displays more clearly the behavior at masses 
where there are too few objects to show on the standard plot. 
At masses below IQ^^Mq, high specific SFRs dominate. 
At M, = 10'°Mo, more galaxies assume values below the 
dominant high-specific SFR "sequence", down to very low 
levels. At = 2 x 1O"M0, close to where the high-specific 
SFR "sequence" terminates, the galaxies with low specific 
SFRs begin to dominate. Above logM, = 11.8 the plot is 
dominated by noise from very few objects with such high 
mass. This figure compares well to Figure 24 in B04. 

7.2. The star-fanning sequence 

We now focus on star formation histories of various classes 
of galaxies. In Figure [17] we show specific SFR against 
the stellar mass for the star-forming, AGN (together with 
SF/AGN composites), and the class without Ha detection. 



For each class the upper panels show nominal, unweighted 
data, while the lower panels are volume-corrected. Dashed 
line has the same meaning as in Figure [15] First we notice 
that the three classes occupy relatively distinct portions of the 
parameter space. This is especially pronounced in the un- 
weighted plots. Thus the three classes appear to have had 
quite different star formation histories. This was to some ex- 
tent indicated in the CMDs (Figure [T]|, but is more striking 
now. 

The star forming class (SF), which forms a blue sequence 
in the CMD (Figure [T]), dominates the high values of the spe- 
cific SFR. The sequence that was obvious in Figure [15] stands 
completely isolated here. In the weighted plot we see an ever- 
increasing density towards the lower mass, which reflects a 
rising low-mass slope of the mass function of late type (blue) 
galaxies. While the typical SFR rises from O.lMQyr"' at the 
low-mass end to lOMoyr"^ at lO'^Mo (consistent with the 
assumption that the more massive galaxies contain more gas 
and therefore have higher SFRs), the specific SFR actually 
declines by a factor of 10. The tightness of the "star form- 
ing" sequence (note that the formal error of the specific SFR 
is only 0.20 dex) represents an important indication that the 
gala xy's mass regulates the overall star forrn ation history (see 
also lGavazzi et alj|1996l: iBoselli et al.ll2001h . It appears rea- 
sonable to assume that in the absence of an event that may 
disturb galaxy's gas reservoir, a galaxy would "naturally" sit 
on this tight sequence. This appears contrary to the sugges- 
tions that the red (non-SF) sequence forms by simple gas ex- 
haustion, since we would then have a much wider range of 
(specific) SFRs at the given mass, for galaxies observed in the 
various stages of gas exhaustion. Note that for masses below 
log = 9.5 we have an additional tail of high specific SFRs 
(best seen in volume-weighted plot), which is indicative of 
starbursts. We will return to the width of the SF sequence in 



7.3. Specific SFRs and stellar masses of galaxies with an 

AGN 

We now turn to galaxies that lie on the AGN branch of 
the BPT diagram. These represent both galaxies classified 
as SF/AGN composites, and "pure" AGN (those lying in the 
region of the BPT diagram where non-SF ionizing source is 
necessary). We have akeady seen that a significant fraction 
of AGN hosts populates the green valley of the UV-to-optical 
CMD (Figure [TJ — the intermediate region between the blue 
and the red sequences. (However, note that in the optical 
CMD they would lie predominantly in the red sequence). A 
possible reason for the presence of these galaxies in the green 
valley is that they are entering the stage in which their star 
formation is being shut off, and the colors are becoming red 
as a result. Martin et al. (2007), using the GALEX data, has 
demonstrated that if the SF of galaxies lying in the green 
valley is indeed being quenched, then their color evolution 
implies a mass flux (transfer from blue to red) that roughly 
matches the mass increase required to populate the red se- 
quence between t he epochs I and z ^ (Bell et al. 2004); 
lFaberetalJl2005h . 

In the middle panels of Figure [17] we present specific SFR 
against mass for AGN and the SF/AGN composites combined. 
The typical error of the specific SFR is 0.35 dex. Galaxies 
with AGN populate the region between the SF galaxies with 
no AGN and the galaxies with no Ha (adjacent panels), in 
a more striking manner than even in the UV-to-optical color- 
magnitude diagram. We note that these AGN are of Type 2 
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Fig. 17. — Dependence of the specific SFR on the stellar mass for different classes of galaxies. Star forming (left panels); galaxies with AGN (middle panels); 
and galaxies without Ha detection (right panels) all occupy distinct regions of the parameter space, indicating different SF histories. SF galaxies form a narrow 
sequence. AGN have intermediate specific SFRs, and are predominantly high mass. Galaxies without Ha, mostly red-sequence galaxies, have low specific SFRs. 
The dashed lines shows the reference SFR of 1 MQyr"' . The lower panel shows values weighted by Vmax- Uneven behavior at low masses is because of small 
number of galaxies (or no galaxies) in those mass bins. 

(Seyfert 2s and LINERs), which are not expected to affect the 
stellar continuum, especia lly not in the integrated galaxy light 
dKauffmann et alj 12003d, and refe rence therein). In galax- 
ies having a luminous Type 2 AGN. lKauffmann et alJ (l2003d) 
limit the AGN continuum contribution to less than a few per- 
cent in r band. Assuming a power-law AGN continuum of 
the form cx i^"'^, this limit corresponds to less than ^ 15% 
contribution in the UV. Another indication that the UV con- 
tribution from an AGN is low in the UV comes from the fact 
that the quality of our SED fitting for galaxies with AGN is 
nearly the same as for the similar galaxies with no AGN. Also, 
iKauffmann et al. (2007) explore GALEX imaging of a more 
nearby sample of AGN hosts and find that their UV emission 
is quite extended, thus unlikely to emanate from the central 
AGN. 

The fact that most AGN have intermediate star formation 
histories, and that this region is not occupied by galaxies un- 
less they have an AGN, is remarkable in the light of recent 
suggestions that the AGN feed back may be regulati ng the star 
formation of massive galaxies (ISpringel et al.ll2065h . Still, we 
would like to know if there is a connection of galaxies host- 
ing AGN with normal star-forming galaxies with no AGN. In 
Figure[T8] we plot SF class (thin solid contours), and separate 
AGN into "pure" AGN (dashed contours) and SF/AGN com- 
posites (bold solid contours). Composites form the lower part 
of the AGN branch in the BPT diagram, and the "pure" AGN 
the upper part. Contours encompass 10, 30, 50 and 70% of 
each class's objects (i.e., their PDF densities), with no Vmax 
weighting. We notice that composites "bridge" the normal 
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Fig. 18. — Specific SFR and mass of star forming, and of galaxies host- 
ing AGN. We plot density contours for normal SF galaxies having no AGN 
(thin solid contours), of SF/AGN composites (bold solid contours), and of 
"pure" AGN (dashed contours). The SF/AGN composite class (lower part of 
the AGN branch in BPT diagram) bridges normal SF galaxies and the AGN 
lying on the upper part of the AGN branch. Contours of unweighted distri- 
bution encompass 10, 30, 50, and 70% of objects, i.e., their composite PDF 
densities. 

SF population and that of "pure" AGN, i.e., they have val- 
ues of both the mass and the specific SFR that are interme- 
diate between the two classes. AGN, in addition, have a tail 
extending to very low values of the specific SFR. Another, 
more physically motivated way of splitting the AGN hosts is 
based on the strength (i.e., the accretion rate) of the AGN it- 
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Fig. 19. — Specific SFR and mass of star forming, and of galaxies host- 
ing AGN of various strengths. We again show SF galaxies with no AGN (thin 
solid contours), galaxies with strong AGN (bold solid contours), and galaxies 
with weak AGN (dashed contours). The strong/weak division of all galaxies 
on the AGN branch of the BPT diagram is based on the [Olll]A5007A lu- 
minosity. Strong AGN appear to form a continuation of the SF sequence, 
having similar SFRs, but larger masses. Weak AGN have similar masses as 
strong AGN but lower (specific) SFRs, with a tail extending to values typical 
of "quiescent" galaxies. Contours of unweighted distribution encompass 10, 
30, 50, and 70% of objects, i.e. their composite PDF densities. 

self, as indicated by its [OIII]A5007A luminosity. We use 
iKauffmann eTal] (l2003cl) dust-corrected [OIII]A5007A lumi- 
nosities, and divide AGN hosts into strong and weak AGN, 
using their demarcation of L[OIII] = IO^Lq. We apply this 
division on all galaxies lying on t he AGN branch of the 
BPT diagram (we refer the reader to IKauffmann et al.ll2003cl 
for a visual presentation of where AGN of various strengths 
lie on the BPT diagram). We show in Figure [19] the re- 
lation between the specific SFR and mass for SF galaxies 
(thin solid contours), strong (bold solid contours) and weak 
(dashed contours) AGN. Apparently, strong AGN lie on the 
massive continuation of the SF sequence, with SFRs as high 
as those of the massive galaxies of the SF class (i.e., the offset 
from the dashed line is similar). Higher masses with similar 
SFRs makes their specific SFRs take more intermediate val- 
ues. Weak AGN fall at lower star formation rates relative to 
the strong AGN, with a tail toward very low star formation 
rates that extend s well into the domain of t he quiescent galax- 
ies (see below). IKauffmann et al.l (l2003cl) . using the optical 
spectroscopy, found that strong AGN have younger popula- 
tions than weak AGN, which we qualitatively confirm, and 
provide a closer connection between strong AGN and normal 
SF galaxies. 

Here, we propose an evolutionary progression for massive 
galaxies: from SF galaxies with no AGN, to strong AGN that 
are more massive yet with similar SFRs, to weak AGN that 
have masses similar to strong AGN but have lower (more sup- 
pressed) SFRs. Thus, the weak AGN would represent a pop- 
ulation with fading star formation. We emphasize that while 
most optical and structural properties of the majority of AGN 
would place them in the optical red sequence, we find that 
they represent a smooth, massive extension of the star forming 
sequence defined by optically blue galaxies. This pic ture is 
supplemented by the results of IKauffmann et an(l2007b . They 
confirm that the majority of blue NUV-r galaxies with high 
velocity dispersion host AGN, and suggest that strong AGN 
require extended gas-rich disks to provide a reservoir of gas 
that is fueled into an AGN. Previously, without the UV data. 



IKauffmann et al.l (l2003d) suggested that many type 2 AGN 
are post-starburst systems. Such a possibility, by itself, allows 
at least two different evolutionary scenarios. One in which a 
normal red-sequence galaxy experiences a burst (say, through 
a minor gas-rich merger) which also fuels an AGN, and an- 
other, in which the AGN is originally concurrent with nor- 
mal SF. Our results, showing a smooth sequence that begins 
with SF galaxies without AGN and extends at its massive end 
to AGN hosts with different levels of star formation, appear 
more compatible with the second s c enario . 

While neither IKau ffmann eldl ( |2007 |) nor our study di- 
rectly proves that AGN are actually responsible for the subse- 
quent quenching of the star formation, they do provide more 
evidence for the connection between the AGN and the fading 
star formation in massive galaxies. Such picture is qualita- 
tively compatible with recent theoretical work based on A^- 
body simulations of galaxy evolution and semi-analytic treat- 
ment of feedback. Specifically, the "ra dio"-mode AGN feed- 
back, proposed by lCroton et al.l (l2006l) . demonstrates that the 
galaxies situated in dark halos above some mass start to sup- 
press the accretion of gas, and at the same time suppress 
further star formation. These are called "radio" AGN be- 
cause they can have black hole accretion rates well below 
Eddington, as is the case in our sample. Also, this mecha- 
nism is efficient in very massive galaxies. From the weighted 
plot of AGN (lower middle panel of Figure [TtI i we find that 
the mass function of AGN peaks at around 4 x lO'^Mo, and 
has an abrupt turn-over a t lower masses. This was noted by 
IKauffmann et al.l (l2003cl) who also demonstrate that the low- 
mass cutoff is not likely to be a selection effect; if a low-mass 
galaxy contained anything but the very weakest AGN it would 
be classified as an AGN in the BPT diagram. The connection 
that GALEX has provided between the SF galaxies and the 
AGN feedback is quite reassuring, and evidently further work 
is required to complete this picture. 

7.4. Specific SFRs and stellar masses of galaxies with no Ha 

emission 

Finally, we concentrate on galaxies for which no Ha emis- 
sion was detected in SDSS spectra. These galaxies have red 
optical colors, which puts them in the "red sequence" cate- 
gory. These galaxies are often referred to as "old, red and 
dead", because of the presumably old populations and the lack 
of current SF. However, in the UV to optical colors (Figure[T]i 
we saw that some of these galaxies have relatively blue colors 
(which indicates that one needs to be careful when using the 
term such as "red sequence", since it apparently depends on 
the color used). In the right-hand panels of Figure[T7]we show 
specific SFRs of galaxies with no Ha. In the unweighted plot 
they occupy a region of the parameter space distinct from the 
star-forming and AGN galaxies, with specific SFRs typically 
lower than for those two classes. The volume-corrected plot 
still puts the bulk of them lower than the other two classes, 
but the overlap with the low specific-SFR tail of the AGN is 
more prominent. Also, the mass extent is now displayed ac- 
curately, with the mass function turning over and flattening at 
lower masses. 

If we expect most red-sequence galaxies to have no star 
formation at all, how do we interpret the (specific) SFRs de- 
rived here? First, we note that the specific SFRs (and also 
just the SFRs) are more poorly constrained than for any other 
class (average formal error in specific SFR is 0.63 dex). This 
should be interpreted in a way that the observed SEDs are 
in many cases compatible with a very wide range of SFRs, 
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mostly very low ones. Let us note that our model SF his- 
tories do not contain truncated models, but only exponen- 
tially declining models (with added starbursts), which can 
never reach exact zero SFR. The dynamic range of SFR in 
models is determined by the range of exponential timescales, 
which in our case have a minimal e-folding time of 1 Gyr 
This produces a lower limit of log(SFR/L) « -14, which 
is roughly equivalent to the lower limit in the specific SFR 
shown in Figure [17] While many galaxies in this class are 
compatible with this lower limit, they do not cluster close 
to it, because in the absence of stronger constraints, higher 
SFRs will also be allowed (i.e., to firmly constrain that the 
SFR is, say, smaller than 0.001 MQyr"' one needs to measure 
a very low amount of UV light with a very high precision). 
However, there are clear cases in which a galaxy classified 
as having no Ha has a relatively well constrained SFR. For 
some of them it is as high as IMoyr"' (crossing the dashed 
line). What is the morphology of galaxies with high (spe- 
cific) SFRs, yet no Ha emission? There are 837 galaxies with 
logM, > 10 and log(SFR/M,) > -11. We inspect SDSS 
images of these galaxies. Most appear as normal early-type 
galaxies (wi th several obvious E-hA types), which was already 
indicated in iRich et aP (l2005l) who find that some morpho- 
logically and spectroscopically early-type galaxies show blue 
UV-to-optical colors. However, one quarter shows some sign 
of structure or a disturbed light profile. This fraction would 
very likely be higher in deeper images or images with higher 
spatial resolution. One quarter of galaxies with structure show 
traces of spiral arms. In most cases these spiral arms are faint, 
and possibly represent cases where the star formation was re- 
cently shut down (leading to the lack of emission lines), yet 
the residual UV light and the spiral structure remain. Such 
galaxies have been recognized in galaxy clusters as "passive 
spirals" jvan den Berghlll976t iPoggianti et al.|[T999h . and are 
believed to represent a transition between a spiral and an SO 
galaxy. In the inspected sample these galaxies are found at all 
redshifts, indicating that the reason Ha is not detected is most 
likely not due to the aperture size being smaller than the size 
of the galaxy from which the integral flux is measured (espe- 
cially for more nearby galaxies). Many of the galaxies with 
structure exhibit rings — from very faint ones to bright blue 
ones. In some of those cases SF could be occurring outside of 
the spectroscopic fiber, leading to non-detection of Ha. 

One can be concerned that our models attribute the UV 
light from old populations (the UV "upturn") to young 
stars, thus implying an ar tificially high star formation 
rate. Bruzu al & Char iot (2003) models do include post- AGE 
phases, but not the extreme blue horizontal branch, which is 
believed (by no means unanimously) t o be responsible fo r 
the UV upturn in early-typ e galaxies jBrown et al.l 120001) . 
However. iDonas et al.l (120071) explore nearby "pure" elliptical 
galaxies and find that both the UV and the optical colors that 
they observe are fully reproduced in the models that we use 
in this work, and that if the UV excess from the actual star- 
formation is present, it is discernible from the combination of 
UV and optical colors. Therefore, our models do appear to ac- 
count for UV light from old populations and therefore some 
of the objects with no Ha emission are quite likely to have 
had actual recent SF. This concl usion is also reached in other 
studies bas ed on GALEX dat a (lYi et al.ll200l iKavirai et all 
|2007; Scha winski et al.ll2007h . where flie blue UV-to-optical 
colors of morphologically selected early-type galaxies are in- 
terpreted as signs of recent star formation. It is quite conceiv- 
able that these galaxies are in a post-starburst phase, so that 



the Ha emission is no longer present, while the traces of star 
formation are still visible in the UV. In evolutionary terms, as 
in the case of AGN, a question remains whether this recent SF 
is the sign of a previously normal SF that is now quenched, or 
is it the result of a small SF episode in an already quiescent 
galaxy. In the latter case, the starburst could be due to a minor 
merger with a gas-rich galaxy. Most likely both processes are 
present, and the forthcoming work will address the relative 
importance of each. 

7.5. Relationships between specific SFR and mass 

One would like to quantify the trends of the specific SFR 
trends with respect to mass in order to facilitate easier com- 
parison of our optically-selected, local-universe sample to 
studies at other redshifts, or samples selected in a different 
manner (keeping in mind that a conversion may be required 
for a different choice of IMF or Hq). 

If one can identify "pure" star forming galaxies (with no 
AGN contribution) using the BPT diagram, then the star form- 
ing sequence is best described with a linear fit: 

log(SFR/M,) = -0.35(logM,-10)-9.83. (11) 

In many instances one does not have all the requisite spec- 
tral lines to make an exact classification, but instead relies on 
color to separate star-forming galaxies. Adopting a color cut 
of NUV - r = 4, which passes in between the blue and the red 
sequence, and taking galaxies on the blue side, we encompass 
all "pure" star forming galaxies, but also many SF/AGN com- 
posites (see Figure [U In that case, the star-forming sequence 
is better described with a piecewise linear fit: 



log(SFR/M,) = 



-0.17(logM, 
-0.53(logM, 



10)- 
10)- 



-9.65, logM, >9.4 
-9.87, logM, <9.4. 

(12) 

This relation differs from Equation [TT] mostly because the 
high-mass end is somewhat lower — a consequence of includ- 
ing SF/AGN composites which have lower specific SFRs 
than pure SF galaxies of the same mass. In both cases we 
have performed least-square fitting to the conditional volume- 
corrected composite PDFs (analogous to the upper panel of 
Figure [16), binned into running modes in the 8.1 < logM* < 
1 1 .7 range for SF galaxies, and 8.0 < log < 1 1 .6 range for 
blue galaxies. 

The trend of the specific SFR with respect to mass for the 
entire sample is best observed in Figure [T6 We find that the 
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parametrization having the form similar to 
luminosity function provides an excellent overall fit to this 
trend. Recently, Feulneretal. (2006) have used Schechter 
functions to describe the upper envelopes of their SFR/M* vs. 

data points, but here we find that it represents a good de- 
scription of the main trend. We use the Schechter function of 
the form: 



SFR /M, = (SFR /M, )o 1 ^' -'°s exp[- 1 0<'°s 'Wo)] ^ 

(13) 

where (SFR/M*)o, Mo, and a describe the normalization, 
characteristic stellar mass, and the low-mass slope, respec- 
tively. Taking modes of the conditional distribution in the 
8.0 < log < 1 1 .9 range, we find the following best-fitting 
parameters: 
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Fig. 20. — Cosmic specific SFR as a function of the limiting mass. We 
define cosmic specific SFR as the the ratio of total SFR to total mass, each 
weighted by Vmax. The figure shows the effect of the limiting mass on the 
estimate of the cosmic specific SFR. SFRs have been corrected to z = 0. 1 by 
applying an evolution correction of the form oc (l+zf"^ ■ 



(SFR/M,)o = 5.96 x 10"" yr"', 

log Mo = 11.03, (14) 

a = -1.35. 

This fit, and the data points that were used to construct it are 
shown in the lower panel of Figure [16] 

In the paradigm of "downsizing", at higher redshifts the SF 
was more active in galaxies with larger mass. This should be 
reflected in the relative flattening of the SF sequence at higher 
redshifts. While most current studies observe an overall shift 
towards higher sp ecific SFRs, th e slope does not appear to 
change much (e.g.. lNoeske et al.] |2007; Papovich et al. 2006). 
Our quantification of the specific SFR should enable more di- 
rect comparison with the local universe. In many cases, sim- 
ple overplotting of Equations [13] and [14] to other sets of data 
would be indicative (unless data are color-selected to include 
only the blue sequence, in which case we recommend using 
Equation[T2]i. 

One often wishes to determine the average specific SFR of 
their sample. The average specific SFR can be defined in sev- 
eral ways, but the most used definition is one in which the 
total SFR in a volume-limited sample is simply divided by 
the total stellar mass. We call this the cosmic specific SFR: 

SFR/M,(Cosmic) = , (15) 

where ° 'SFR is the SFR corrected to z = 0.1 (see ®. From 
our sample we find SFR /M, (Cosmic) = 6.8 x 10"" yr"', 
which is in very good agreement with the equivalent measure 
in B04. How much is the estimate of the cosmic specific SFR 
dependent on the limiting mass of the sample? In Figure l20l 
we plot cosmic specific SFR as it would have been determined 
had we had a mass limit (and assuming the same magnitude 
limits). Right-handy-axis shows a fraction of the "actual" (no 
mass limit) cosmic SFR/M, at the given mass limit. 

Finally, we ask with what precision can one use the rela- 
tions to "predict" specific SFR from stellar mass? For blue 
galaxies, we find the average scatter around the mean se- 
quence to be 0.5 dex. This result is only slightly affected after 
taking out in quadrature the formal measurement error of 0.2 
dex. Therefore, the sequence has an intrinsic width, and is 



TABLE 2 

Comparison OF SFR densities at z = 0.1 in this study and in 

B04, DERIVED USING COMPATIBLE METHODOLOGY 



Class 


PSFR 


% 


PSFR 


% 


fractional 




This paper 


of total 


B04 


of total 


change 


All 


1.953 


100.0% 


1.980 


100.0% 


0.99 


SF 


1.395 


71.4% 


1.013 


51.2% 


1.38 


SF (low S/N) 


0.296 


15.2% 


0.455 


23.0% 


0.65 


Comp 


0.147 


7.5% 


0.207 


10.5% 


0.71 


AGN 


0.076 


3.9% 


0.084 


4.2% 


0.90 


No Ha 


0.038 


2.0% 


0.221 


11.2% 


0.18 



Note. — For this comparison we adopt /3 = 3 SFR density evolution. 
SFR density is given in units of 10""/i7oMQyr"'Mpc"^ 



likely to depend on other factors besides the mass. The width 
itself is a function of mass. It declines at the rate of -0.11 
dex"', i.e., the scatter of the specific SFR at logM* = 8.0 is 
65% larger than at log = 10.5. This is an indication that the 
star formation becomes less stochastic as the mass increases, 
which is what one expects considering that massive galax- 
ies have a larger number of star-forming regions than dwarf 
galaxies. 

8. GLOBAL STAR FORMATION DENSITY 

In this section we use our SFRs to determine the global star 
formation density at z = 0.1, the mean redshift of our sam- 
ple. In order to take into account that the sample has a spread 
of redshifts, and therefore probes a range of cosmic star for- 
mation history, in the remainder of the section each galaxy's 
SFR estimate is "corrected " to the mean redshift of the sur- 
vey, where the uncertainties in the correction are minimized. 
To make the correction, one assumes a cosmic star formation 
density evolution of the form (1 +1}^^, where B04 take /3 = 3. 
Therefore, the corrected SFR for a galaxy lying at redshift z 
is " 'SFR = [1.1/(1 +z)]^SFR. 

Since we will be comparing our results to those obtained by 
B04, we first describe their methodology and apply it to our 
data. B04 obtain star formation density (psfr), and its for- 
mal (random) error, by Monte Carlo method applied to SFR 
PDFs in the following way. First, 100 different bootstrap sam- 
ples are constructed from the original sample, using drawings 
with replacements. Bootstrap samples have the same num- 
ber of objects as the original sample. A random number is 
then used to select a value from the cumulative PDFs for each 
galaxy in the bootstrap sample. The random sampling of the 
PDF is repeated 30 times. Each SFR PDF is weighted by 
1 /Vmax and summed to get psFR- From the ensemble of such 
samplings the mode and the 16% and 84% confidence lev- 
els are determined for the entire sample, and for each galaxy 
class. The results of this exercise are presented in B04 Table 
3. The random error of the total psFR is found to be around 
0.018 in units of lO"^/z7oM0yr"'Mpc""'. For our sample, we 
perform only the multiple probing of each galaxy's cumula- 
tive PDF, without bootstrapping the sample. Performing 100 
Monte Carlo summations, we find a random error of the to- 
tal psFR to be 0.017, in the same units. Therefore, the two 
studies have comparable formal random errors, which is not 
surprising given the fact that the formal errors of the SFRs are 
comparable. In any case, these random errors are smaller than 
some systematic uncertainties that we will explore later 

How do the actual values of SFR density compare? First, 
we notice that the B04 method is essentially equivalent to the 
summation of the medians of the PDFs since the average value 
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TABLE 3 
SFR DENSITIES AT z = 0.1 
DETERMINED IN THIS STUDY 



Class 


PSFR % 


of total 


All 


1.828 


100.0% 


SF 


1.328 


72.6% 


SF (low S/N) 


0.277 


15.2% 


Comp 


0.136 


7.5% 


AGN 


0.067 


3.6% 


No Ha 


0.021 


1.1% 


NOTE. — We adopt 13 = 2.5 
SFR density evolution. SFR 
density is given in units of 
lO"^/i7oM0yr'Mpc"^. Random 
error of SFR density for 'AH' galax- 
ies is 0.017 in same units. See text 



for the discussion of systematic 
errors. 

of a random number that probes the cumulative PDF is 0.5. In 
Table |2] we compare the values obtained from the summation 
of the medians of our SFR PDFs, to the average of the 16 
and 84 percentiles of B04 values (which should be close to 
the median). Our estimate includes the correction for galax- 
ies excluded from the sample because of the poor SED fits 
04.21 1. These corrections are 1.6-2.8%, depending on the 
class. Looking at the first row, we see that our estimate of 
the total SFR density is very close to that of B04. However, 
this is a mere coincidence when we look at how the densities 
compare by class. Our SFR density for the SF class is 38% 
higher; however, for all other classes the estimate is lower — 
between 10% less for AGN to 5.7 times less for galaxies with- 
out Ha. The causes for these differences can be traced from 
the analysis of SFRs in |5] Consequently, the contribution of 
each class to the total SFR density has changed as well. SF 
class now contributes some 70%, while it was around 50% in 
B04. This comes at the expense of other classes (except for 
the AGN), and most significantly from much lower contribu- 
tion of galaxies without Ha (2%, vs. 11% in B04). 

We next describe the calculation of our SFR density, which 
somewhat differs from B04. Direct summation that was used 
to find SFR density in B04 assumes that the losses from the 
incompleteness of the sample are insignificant. B04 estimates 
them to be between 1 and 2 percent. We fit a Schechter fit 
to the number density of galaxies of various classes and inte- 
grate the luminosity density down to our completeness limit 
of Mr = -14.75, and to a limit which is 10 magnitudes fainter. 
We compare the two total luminosity densities and find that 
1.1% of luminosity density is below our limit for SF (and 
low-S/N SF) galaxies. For other classes, the loss is practi- 
cally zero. Therefore, for our SFR density estimate we per- 
form the direct summation by taking galaxies M,- < -14.75, 
and then correcting the SFR density. This result differs in- 
significantly from one in which we would sum all galaxies 
with no magnitude limit. Our SFR density estimates are ob- 
tained by summing the averages of the SFR estimates. We 
prefer summing averages over medians because it is equiv- 
alent to summing an entire PDF. Averages produce a result 
that is 5% higher than the medians. Also, we adopt the expo- 
ne nt of the SFR dens ity evolution /3 = 2.5 ±0.7, determined 
by lSchiminovicheta l. (2005) from GALEX-VYDS measure- 
ments. 

The results, and the breakdown by class, are presented in 
Table [3] and shown in Figure 1211 These are our nominal SFR 




SF Sl<'(IowS/N) Comp AGK Ko Ha 

Fig. 2 1 . — Total star formation rate densities for different classes of galax- 
ies. 88% of total stai' formation occurs in "normal" SF galaxies (including 
low-S/N SF class), and another 1 1% in galaxies that host an AGN. Up to 1% 
occurs in galaxies that are quiescent based on the optical spectroscopy (no 
Ho class). Estimates are given for the mean redshift of the sample of z = 0. 1, 
by applying an evolution correction of the form cx: (l+z)^'^. 

densities. Galaxies for which no classification was possible 
(see ^3.2| l are included in the SF class. They constitute 0.8% 
of the total SFR densuty. Altogether, 88% of star formation 
occurs in "normal" SF galaxies, and another 1 1 % in galax- 
ies that host an AGN. Up to 1% occurs in galaxies that are 
quiescent based on the optical spectroscopy. 

The uncertainty in the evolution exponent leads to the sys- 
tematic uncertainty in the total SFR density of 2%. By far the 
largest uncertainty stems from the systematic uncertainties in 
the determination of the dust attenuation, and hence the SFR 
of SF galaxies ( ^5.4l i. We estimate the degree of this error by 
simply replacing the SFRs of SF galaxies by those obtained 
by B04. This error is preferentially in the direction that the 
"Ha" estimate becomes 1 1 % higher for the SF class, lead- 
ing to the total density that is 8% higher We also take into 
account that our SFR estimates are averaged over 10*^ yr. Us- 
ing 10^ yr instead, we obtain a SFR density that is 1% lower, 
which we will treat as the systematic error. Finally, another 
small source of systematic error comes from the error in the 
estimate of survey area of 0.3%. Adding up in quadrature ran- 
dom and systematic errors in the positive and in the negative 
directions, we arrive at our total, dust-corrected SFR density 
estimate at z = 0.1: 

PSFR = 1.828![!:i^^ X 10-2/27oMoyr-■Mpc-^ (16) 

which assumes the Chabrier IMF. For a compila tion of recent 
measurements of the local psFR please refer to ' Hanish et alj 
(2006). Our psFR agrees very well (< la) with four of the six 
estimates listed in their Table 1, but has significantly smaller 
error bars. Note that except for B04, all other estimates are 
based on samples with fewer than 300 galaxies, and many ap- 
ply fixed dust correction for all galaxies. We can find the aver- 
age FUV attenuation by "de-correcting" our SFRs by adding 
our FUV attenuation estimates to each galaxy. Such uncor- 
rected SFR density is smaller by a factor of 6.34. In other 
words, the average FUV attenuation for the volume-complete 
sample is Afuv = 2.01 mag. 

We next explore the distribution of global SFR density 
across the galaxies having different properties. In Figure 1221 
we show the distribution of SFR density against the rest-frame 
NUV - r color. The distribution sharply peaks at ^\NUV-r) = 
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Fig. 22. — Distribution of tlie star formation density (at z = 0.1) as a 
function of color. We choose UV-to-optical color which shows a distinct 
bimodahty in the number density of galaxies. The SFR density distribution is 
unimodal, and peaks sharply. There is a tail towai'ds the red galaxies, which 
is in part due to the intrinsic dust reddening. A large fraction (94%) of the 
star formation density arises from galaxies blueward of NW -r = 4 — our 
nominal division between the blue and red galaxies. First bins includes all 
objects without NUV detection. Objects below our absolute magnitude com- 
pleteness limit (~ 1% of SFR density) are not included. The units of SFR 
density distribution are 10"-/i7oMQyr"'Mpc"'mag"' . 
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Fig. 23. — Distribution of the star formation density (at z = 0.1) in galaxies 
with different specific SFRs. Specific SFR is an indicator of a galaxy's star 
formation history. The distribution shows a "bump" at the high end due to 
the population of bursty, low-mass galaxies. There is a steep decline towards 
the galaxies with low specific SFR. The units of SFR density distribution are 
10"-/!7oMQyr"'Mpc"^^dex"' . 

2.2. While the number density distribution is strikingly bi- 
modal in this color, the SFR distribution is very even. The 
blue edge is more abrupt then the red side, which shows a tail 
extending to red-sequence colors, in part due to the intrinsic 
dust reddening. We have previously used an NUV -r = 4 cut, 
which goes across the green valley, to conveniently divide the 
blue and the red sequences. From Figure |22] we see that the 
large fraction (94%) of star formation occurs blue of this cut. 
The first bin sums SFR density in galaxies lacking > 3(7 mea- 
surement in NUV (2% of total). 

UV to optical color is easy to measure, but gives only a very 
rough idea of a galaxy's SF history. Therefore, in Figure |2?] 
we plot the SF density distribution against the specific SFR. 
The picture is somewhat different than from the distribution 
against the color First, the high specific-SFR end does not 



Fig. 24. — Distribution of the star formation density (at z = 0.1) in galaxies 
of different stellar mass. The fractional star formation density increases with 
galaxy's stellar mass and reaches a maximum at ai 4 X IO'^Mq, which is 
close to the "transitional" mass of [Kauffmann et al. l2003b). Our sample 
is incomplete at log M» < 7.5. The units of SFR density distribution are 
lO-2/!7oM0yr-' Mpc-^dex"' . 
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Fig. 25. — Fraction of the total star formation density (at z = 0.1) as a 
function of hmiting stellar mass. The summing is performed from the high- 
mass end. Therefore, at each mass, the plot shows how much SF is missed by 
not including less massive galaxies. Our sample is incomplete at log M* < 
7.5, where we estimate 1% of SF takes place. 

fall as evenly as did the blue NUV - r color end; instead, we 
see an excess of galaxies with log (SFR /M*) « -8.3. These 
are the bursty, low-mass galaxies mentioned in ^7.21 The rea- 
son that NUV - r color does not distinguish this population is 
because it "saturates" for galaxies with log(SFR/M,j) sa -9, 
so that all galaxies with higher specific SFR will have sim- 
ilar NUV - r colors. FUV (included in the specific SFR) is 
needed to break the degeneracy. On the side of the low spe- 
cific SFR (equivalent to redder A^f/V-r colors), the decline is 
more steep than for the colors. The difference arises from the 
fact that the specific SFR, unlike NUV - r color, does take out 
the effects of the intrinsic dust attenuation. 

Finally, in Figure |24] we show the distribution of SFR den- 
sity across the galaxies with different stellar masses. The frac- 
tion of SF occurring in galaxies with increasing mass rises al- 
most linearly. The peak of the SF today occurs in galaxies 
having a mass of w 4 x lO^^Mp^. This mass is close to the 
"transitional" mass of iKauffmann et al] ^2003b), which is ef- 
fectively a mass at which blue and red mass functions cross. 
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i.e., the mass above which the red galaxies become more nu- 
merous. In the paradigm of downsizing, we would expect this 
distribution to move towards the higher masses at higher red- 
shifts. This, however, is often difficult to establish since the 
completeness limits of many high-redshift surveys are uncom- 
fortably close to the peak of the distribution whose position 
one tries to establish. In Figure |25] we show the cumulative 
SFR density as a function of stellar mass, summing up from 
the highest mass galaxies. This figure shows how much SF 
is missed (in the local universe) by being limited to a given 
stellar mass (and assuming the same magnitude limits). Thus, 
a limit of lO'^M© captures only 6% of the total SFR density. 
In order to reach 50%, one needs to probe galaxies more mas- 
sive than 1O'°M0. Today, 1/2 of the total SFR density comes 
from galaxies in the mass range 9.3 < log < 10.6. 

9. CONCLUSIONS 

We have considered the measurement of the star formation 
rate using a sample of 50,000 optically selected galaxies with 
GALEX and SDSS photometry and SDSS redshifts. We have 
given a detailed description of our approach to SED fitting, 
including the use of a large library of population models that 
include a physically motivated treatment of dust attenuation. 
These give dust-corrected SFRs ("UV" SFRs) of galaxies in 
the local universe. We show that there are no inherent prob- 
lems that would prevent UV from serving as a reliable SFR 
indicator, despite past concerns over population synthesis and 
dust attenuation models. We recite some of the highlights 
from our study. 

• 25 min imaging from GALEX in two UV bands can be 
used to derive SFRs with a formal accuracy of 0.2 dex 
for r < 1 8 mag star-forming galaxies with known red- 
shifts. 

• For galaxies for which the dominant ionizing source 
are young stars, (as determined using 3" fiber mea- 
surements and the BPT diagram) the SFRs obtained 
through SED fitting ("UV" SFRs) agree exceptionally 
well with dust-corrected SFRs obtained from modeling 
of the nebular lines — essentially Ha, corrected using 
a Balmer decrement ("Ha" SFRs), as long as the same 
dust model is applied to both methods. The comparison 
presented here is at a scale that is two to three orders of 
magnitude larger than the comparison of these SF indi- 
cators done in the past. 

• The lCharlot & Falll ( l2000 l) dust atte nuation model that 
we ap ply to "UV" SFRs and that iBrinchmann et al.l 
(12004 ') apply to "Ha" SFRs succeeds in producing UV 
and Ha dust attenuation estimates that lead to compat- 
ible SFRs. 

• Minor differences between the "UV" and "Ha" SFRs 
exist at the extremes of galaxy stellar mass distribu- 
tions (equivalent to the extremes of the dust attenua- 
tions). They arise either from somewhat weaker con- 
straints from the SED-derived dust attenuation, or from 
the mass-dependent deviations of the assumed oc A"°^ 
extinction law in the UV regime. 

• If the systematics in the attenuation estimates are ac- 
counted for, the "UV" and "Ha" SFRs agree within 
10% across the entire range of galaxy masses. 



• We calibrate the relations that allow SFRs to be deter- 
mined directly from K-corrected FUV and NUV mag- 
nitudes. They involve an FUV attenuation estimate 
from the UV color (equivalent to the UV slope) that we 
calibrated for normal sta r-forming gala xies, and find it 
to be less steep than the iMeurer et alJ (1999) relation. 
To obtain SFRs we find that one requires a conversion 
factor between t he FUV lu minosity and SFR that is 
30% lower than iKennicuttI ([1998) value. Such "sim- 
ple" UV SFRs compare well with those obtained from 
the full SED fitting in the 0. l-lOMoyr"' range, i.e., for 
normal star-forming galaxies. 

• For some 50% of the optically selected galaxies from 
SDSS, the emission lines are either too weak (red se- 
quence galaxies), or they are contaminated by narrow- 
line AGN, which precludes obtai ning SFRs from Ha 
directly. In those cases Brinchmann et al.l (120041) use 
calibrations between SFR and optical colors from SF 
galaxies, and apply them to other classes of galaxies. 
We find that such calibrations are not reliable and lead 
to large overestimation of SFRs. UV, on the other hand, 
provides more secure SFR estimates. 

• We present stellar-mass dependent star formation his- 
tories for various classes of optically selected SDSS 
galaxies, from the starbursting dwarf galaxies with spe- 
cific SFRs (SFR/M*) of lO'^yr"' to essentially dead, 
red-sequence galaxies with 6 orders of magnitude lower 
specific SFR. We also present a volume-corrected ver- 
sion of these relationships which give a better account 
of the SF at low masses. We quantify these relations in 
order to facilitate comparison with similar results ob- 
tained at higher redshifts. 

• Star forming galaxies with no AGN form a very well 
defined linear star-forming sequence in specific SFR 
vs. mass plot, with the intrinsic 1 <t width of 0.5 dex 
in log SFR/M*. The sequence is 65% wider at the 
low mass end (^ IO'^Mq) than at high masses, consis- 
tent with a more stochastic (bursty) SF history for gas- 
rich dwarfs. The star-forming sequence terminates at 
3 X 1O"M0. 

• Galaxies whose ionizing flux (as measured in 3" fiber) 
comes primarily from young stars (star formation), 
those for which there is a narrow-line AGN component 
(mostly Seyfert 2s and LINERs), and those which are 
spectroscopically quiescent (no detection of Ha lines), 
all occupy distinct positions in the specific SFR vs. stel- 
lar mass diagram, with the AGN hosts lying in between 
SF galaxies with no AGN, and the quiescent galaxies. 
Star forming galaxies with no AGN do not occupy the 
intermediate regions of the specific SFR — those are oc- 
cupied exclusively by galaxies with AGN. 

• If we take all galaxies forming the AGN branch in 
the BPT diagram (again from fiber measurements) and 
divide them by their AGN strength (as indicated by 
[OIII]A5007A luminosity), we find that strong AGN 
represent the massive continuation of the SF sequence 
defined by galaxies with no AGN. Strong AGN have 
similar SFRs as the most massive galaxies with no SF, 
but are more massive, driving their specific SFRs to in- 
termediate values. Weak AGN have similar masses as 
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strong AGN, but fall at lower star formation rates rel- 
ative to strong AGN, with a tail toward very low star 
formation rates that extends well into the domain of the 
optically red galaxies with no Ha. While AGN appear 
to be a massive continuation of the blue, star-forming 
sequence, their optical colors would place most of them 
onto the red sequence. We suggest an evolutionary se- 
quence from massive star forming to quiescent galax- 
ies via strong and weak AGN that bridge them. This 
scenario supports a picture of AGN feedback suppress- 
ing star formation in galaxies situated in very mas- 
sive dark matter ha los, the so-called radio mode AGN 
dCroton et al.ll2006l) . Further study is needed. 

• While mostly inactive, some galaxies with no de- 
tectable Ha in their fiber measurements exhibit SFRs as 
high as ^ IMQyr"'. The optical colors of these galaxies 
are predominantly red. We find that many (1/4) of such 
galaxies display disturbed light profiles, sometimes in 
the form of fading spiral arms with no current SF de- 
tectable in Ha. In many cases we see ring-like struc- 
tures. Some of these are likely post-starburst systems. 
We find that up to 1 % of U V SFR density comes from 
galaxies with no Ha lines. 

• We determine total SFR density in the local universe 
iz = 0.1) of 1.828!!]:(i^[^ X 10-2/i7oMoyr-'Mpc-^ (for 
Chabrier IMF), which is in good agreement with most 
recent estimates, yet has a smaller error (which includes 
random and most systematic errors). 

• Large fraction (95%) of SF takes place in blue sequence 
galaxies (NUV - r < 4). SFR density increases in 
galaxies of increasing mass, and peaks in galaxies with 
m asses 4 x IQ^^^Mp, — clo se to the "transition" mass 
of iKauffmann et al] (f2003bi) (i.e., the mass at which the 
mass functions of red (old) and blue galaxies intersect). 
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Galaxies in the mass range 9.3 < logM, < 10.6 ac- 
count for 1/2 of the total SFR density. 

We demonstrate the ability of GALEX UV observations of 
the galaxies in the local universe to constrain SFRs of a sta- 
tistically large number of galaxies with a wide dynamic range 
of SF activity. Besides characterizing the galaxies in the lo- 
cal universe, this study provides a better reference point for 
high-redshift investigations. This and the forthcoming stud- 
ies based on GALEX will continue to provide a better under- 
standing of the galaxy evolution, specifically the origin of the 
galaxy bimodality, the inter-connection between the SF and 
AGN activity, and the general star formation history of galax- 
ies spanning a range from gas-rich dwarfs to massive ellipti- 
cals. 
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